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attrition must be based upon sound mathematical formulations and parameter esti- 
mations. With an inherent lack of representative combat data for modern warfare 
scenarios, one effective method for determining the required parameter estimates is 
to thoroughly analyze the output from a stochastically based high-resolution combat 
model. It is this development of attrition parameters process, which so profoundly 
influences the validity of aggregated simulations, that lacks any comprehensive doc- 
umentation or mathematical justification within the modeling community. By ex- 
amining the development and validity of these processes for parameter estimation, 
valid attrition calibration formulae can be determined and used within force attri- 
tion algorithms in order to more precisely and justifiably model aggregated combat 
operations. The establishment of a user-friendly test bed for examining this attri- 
tion rate development process will play a major role in solidifying the understanding, 
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i INTRODUCTION 


Combat modeling and simulation are important analysis tools available to mil- 
itary leaders and their supporting organizations which provide valuable insight into 
the study of military operations across the spectrum of conflict. Combat models 
are employed frequently to conduct investigative studies ranging from peacekeeping 
operations, which are relatively narrow in scope, to conventional maneuver warfare 
executed at the nation-state level. Commanders and their supporting planning staffs 
regularly seek the resulting output from the execution of combat models and simu- 
lations. Prior to the employment of man and machine into a hostile environment, 
tactical commanders have the capability to develop and assess their maneuver plans 
by use of combat models and simulations. Additionally, combat models and simu- 
lations are employed as training tools to rehearse command and control actions and 
operational procedures. While models and simulations are not in themselves predic- 
tors of the outcomes to military operations or plans, their ability to provide samples 
about possible outcomes is welcomed as a preliminary step prior to placing soldiers 
and equipment in harm’s way. As such, the designers of these models constantly 
attempt to refine the representation of warfare to ensure the model is as accurate 
and believable as possible. It is important to note that many of the interactions and 
human factors which are inherent to a military force are difficult to replicate within 
a combat model or simulation. This fact must be recognized by those employing 
combat models and simulations for mission analysis and maneuver warfare planning. 

One of the most important capabilities of combat models and simulations is 
their ability to calculate the attrition of personnel and systems of the engaged military 
forces. The prediction, analysis, and management of the attrition to personnel and 
equipment in combat are highly regarded as important tasks by military leaders and 
their staffs. And, as a result, the ability of combat models and simulations to capture 


attrition data is often a measurement of the model’s usefulness. What follows is a 


discussion of the development of attrition measurements within the context of combat 


models and simulations. 


A. BACKGROUND 


Several important principles and concepts underlying the general class descrip- 
tion of combat models and simulations are important to this development of attrition 
measurement. Combat models and simulations are generally subdivided into two 
categories, namely, high resolution and aggregated. 

High resolution models are generally associated with smaller force compo- 
sitions where each combatant and its associated weapon or maneuver system are 
explicitly represented in the simulation. The complicated combat processes are de- 
composed into highly detailed sequences of events and activities which are carried 
out explicitly within the simulation. This detailed sequence of events and activities 
lends itself to simulation time steps which are very short. It is not uncommon for 
high resolution models to represent time explicitly to the second, updating 1mpor- 
tant information almost continuously. The goal of the high resolution model is to 
accurately model each combat process so that the simulation results are traceable 
using mathematical formulas and constructs. The characteristics associated with the 
entities represented in the simulation are related to their engineering specifications, 
physical performance data, and specified behavioral assumptions. The existence of 
this information link from the engineering level through the behavioral engagement 
assumptions is the single most important quality of a high resolution model. The 
capability of representing weapon systems and equipment at the engineering level 
enables an analyst to observe the possible effects of alternative system specifications 
on a pre-defined simulation. Individual system contributions to the overall simulation 
are effectively analyzed in high resolution models. Another feature associated with 
high resolution models is based upon their reliance on stochastic processes. Nearly 


all of the simulation algorithms above the engineering specification level are proba- 


bilistic in nature. For example, consider a firing sequence. The firing system goes 
through several sequential actions prior to firing its weapon. These include target 
acquisition, target selection, firing, impact, and damage assessment. As these events 
are not of a fixed length in actual combat, their representation in simulation relies 
upon the use of highly probabilistic random variables. Calculations of this nature 
can lead to variance among independent trials of a given combat simulation. In order 
to obtain a reliable result with respect to the variance between stochastic calcula- 
tions, the modeler must execute multiple runs of the same simulation to obtain the 
“averaged” result. One of the most significant limitations of high resolution models 
is their sheer size and complexity. The computer codes supporting high resolution 
models are highly complex and often quite large, stretching hardware and software to 
their limit. These computer codes are expensive to develop, run, and maintain. The 
complexity associated with representing each individual entity of the combat scenario 
at the engineering specification level causes another significant limitation of high res- 
olution models, namely, force size. High resolution models are limited to relatively 
small forces. In ground maneuver terms, this limitation is normally associated with 
an upper bound of roughly 1000 combatants [Ref. 1] or on the order of battalion to 
brigade levels of combat. 

Aggregated models, as the word implies, maintain varying levels of aggregation 
of the represented combatants. For example, an aggregated model may be capable 
of representing a group of riflemen with an associated consolidated firepower capa- 
bility rather than the individual soldiers. The model then portrays the aggregated 
capabilities of the group of riflemen represented graphically as a single entity. There- 
fore, an entity representing a platoon of 30 riflemen would maintain the firepower 
and maneuver capabilities associated with the 30 individual weapons systems carried 
by the 30 riflemen. The extent of aggregation within combat models varies based 
upon the level of combat being modeled and the intended use of the model. The 


aggregation process is natural within the confines of military combat models because 


of the hierarchical command structure inherent to maneuver forces. For example, 
the natural aggregation process associated with a generic ground maneuver organiza- 
tion is individual combatant, squad, platoon, company, battalion, brigade, division, 
corps, army, and theater. Depending upon the use of the model in question, one 
of these natural groupings (normally no lower than brigade) would mark the lowest 
level of representation within the simulation. Individual combatants and groups of 
combatants which are smaller than the basic size are not explicitly represented and 
the attributes associated with these combatants are lost to the simulation. ‘Two im- 
portant and distinctive aggregation patterns associated with models of this type are 
homogeneous and heterogeneous aggregation, and, the distinction between them is 
particularly important for attrition calculations. A homogeneous aggregation is one 
where the combat power of a unit with multiple and distinct combatants is combined 
into a single measure (or perhaps one for ground combat systems and a second for 
aircraft systems). Attrition computations are then based upon the relative power of 
the aggregated forces in battle subject to their combat power indices. Contrastingly, 
heterogeneous aggregation maintains a count of the number of surviving weapons sys- 
tems of each distinct type. This method of aggregation allows the model to represent 
the different levels of effectiveness of a particular weapon system against its possible 
targets. This form of aggregation permits more accurate attrition modeling than the 
homogeneous model resulting in less information loss. As a result, the trend of most 
modern aggregated models is to follow a heterogeneous aggregation process. 

As the size of the modeled entity increases in aggregated models, the preci- 
sion to which they can be located on the simulated battlefield may be decreased. A 
modeled combatant may maintain a center location and a frontal width where the 
individual entities making up the aggregated unit are assumed to be somewhere in the 
unit region. Information about the individual entities within the aggregated unit may 
be kept in the form of averages, but any specificity about individual locations is lost. 


Thus, the first significant consequence of aggregation is the loss of individual char- 


acteristics within aggregated units resulting in the use of their averaged properties 
where any inherent fluctuations have been smoothed. Analogous to the smoothing 
of individual variations is the handling of stochastic processes. Since aggregate mod- 
els do not maintain a reference to what each individual combatant is doing at any 
given time, the complex probabilistic events found within high resolution models are 
replaced by averaged behaviors. The firing process described earlier (target acquisi- 
tion, target selection, firing, impact, and damage assessment) within the context of 
high resolution models may be replaced in an aggregated model by computing the 
rate at which an average tank kills enemy tanks. This derived attrition rate implic- 
itly defines each of the probabilistic processes described above with the stochastic 
variability suppressed. This fact leads to a descriptor generally associated with ag- 
gregated models, namely, deterministic. With the engineering specification level of 
detail and the probabilistic processes being suppressed within the aggregated model, 
many assumptions and scenario details get absorbed by the attrition rate number, 
which in itself is a highly variable entity. 

Lastly, aggregated models lose the information associated with event sequenc- 
ing since they do not keep track of individual actions. Thus, the precise times of 
critical events, such as target kills, are not available from aggregated model output. 
This results in the process averages being processed over relatively long periods of 
time (ranging from a few minutes to a 24 hour period) where the computation of total 
casualties occurs at the end of the time period. Again, careful consideration must be 
given to the attrition rate formulation in that there are at least two sides involved in 
the exchange of fires. Not all combatants on either side will survive to continue firing 
for the duration of the time period. For example, consider the following example and 


resulting conclusions: 


Extreme care must be taken to make the underlying assumptions clear 
in computing the rate data for such models. For example, a daily attrition rate 
for aircraft stationed at an airbase which is under attack may be inconsistent 
with the day’s plan for those aircraft (they may be flying missions and thus 


not on the ground when their base is attacked). A simulation which assesses 
attrition on an aggregated daily cycle cannot know whether the airbase attack 
coincided with aircraft on the ground, so an assumption must be made and 
figured into the attrition rate data. 


Overall, then, the outstanding feature of aggregation in combat models 
is information loss. As a result, it is necessary to model combat processes using 
average behavior rather than individual behavior. Although values for the 
average process rates are generally computed outside the combat model, the 
complex nature of these rates demands that we devote considerable effort...to 
explaining their derivation. [Ref. 1] 


An obvious question with regards to the given strengths and limitations of the 
described combat modeling and simulation approaches is why high resolution models 
are not extended to larger force sizes thus maintaining the admirable qualities of the 
high resolution definition. Unfortunately, the difficulty in making this extension lies 
within the limitations of computer hardware and software and the expansive execu- 
tion time which would be associated with such large scale high resolution models. 
The sheer number and complexity of the individual combatants and their associated 
characteristics render the use of high resolution combat models for all combat simu- 
lations impractical and, in some cases, impossible. Additionally, the complex nature 
of command, control, and communications (C*) between echeloned organizations is 
extremely difficult to model. The synergism and inherent lack of measurability of 
these traits would require many assumptions on the part of the modeler possibly re- 
ducing the degree of model believability. Thus, the resulting segregation of combat 
models into the high resolution and aggregated categories is a matter of necessity. 
It is important to note, however, that the two categories may not be entirely inde- 
pendent during their execution due to the aggregated models’ possible reliance upon 
the outcomes of high resolution battle simulations for the derivation of attrition rate 
parameters. While there are several currently employed methods for calculating force 
attrition within aggregated models (to be presented in Chapter III), it remains true 


that the attrition of personnel, equipment, and supplies resulting from an engagement 


must, at some point, be calculated within an aggregated simulation. The compromise 
which is of specific interest to this research is the algorithmic application of mathe- 
matical constructs to determine aggregated force attrition based upon high resolution 
results. For example, when an aggregated combat simulation reaches a point where 
two opposing forces engage in battle, the model may implicitly rely upon high res- 
olution simulation results generated under similar environmental, force composition, 
and tactical conditions. The end result of the battle executed within the aggregated 
simulation would be new force levels for each of the engaged combatants derived pri- 
marily from a deterministic replay of a representative battle using the attrition rate 
coefficients derived from the high resolution simulation output. The intermediate 
process of estimating the attrition rates and replaying the battle mathematically is 
the conceptual link between high resolution and aggregated models. With the line di- 
viding high resolution and aggregated models clear in definition only, both categories 
have an associated collection of applications which, when combined, adequately meet 


the needs of combat modeling and simulation. 


B. APPLICABILITY TO CURRENT PROCEDURES 

In the United States Army and, moreover, the joint warfare community, one 
of the most frequently used high resolution models which provides the data necessary 
to generate the attrition rates for use in aggregated models is called the Combat 
Sample Generator (COSAGE) model. The COSAGE model is a (fairly) high reso- 
lution, stochastic model of division level ground combat which maintains a combat- 
ant resolution to friendly platoons and enemy companies. Within each maneuver 
unit, a heterogeneous list of individual weapon systems is arranged in combat forma- 
tions, interactions between opposing weapon systems are computed, and individual 
weapons system survival is governed by the laws of probability. With respect to 
aircraft, close air support of the ground battle by tactical aircraft is resolved down 


to a detailed simulation of the flight of individual aircraft. At a point in the exe- 


cution of an aggregated model where two opposing ground forces are positioned to 
engage in combat, the representative battle simulation is carried out as a mathe- 
matical algorithm based primarily upon the results of a high resolution simulation 
(COSAGE). The killer/victim results and the associated time series of individual 
attrition are analyzed and mathematically manipulated resulting in the development 
of attrition rate estimates for each individual weapon systems of each engaged force. 
This process is referred to as the attrition calibration (ATCAL) Phase I process and 
will be discussed in detail in Chapter IV. These derived attrition rates are then 
applied to a deterministic algorithm which, when executed or “replayed” with the 
force levels present in the aggregated model, calculates resulting force levels (by 
system). These updated quantities are then passed back to the aggregated model 
as it continues through its remaining processes. This secondary “replay” process 
is commonly referred to as ATCAL Phase II. Currently, several theater level mod- 
els have the capacity to employ ATCAL Phases I and IJ]. These include Tactical 
Warfare (TACWAR), Concepts Evaluation Model (CEM), and Integrated Theater 
Engagement Model (ITEM). Whether a particular theater level aggregated model 
employs the ATCAL algorithm as its method for calculating attrition is largely de- 
pendent upon the respective developers’ and users’ understanding of ATCAL’s theo- 
retical development and mathematical principles, which, in themselves, remain to be 


areas of active research. See {Ref. 2]. 


C. PROBLEM STATEMENT 


In order to support current and future investigations into the theoretical un- 
derpinnings and mathematical stability of the defining algorithms within ATCAL 
Phases I and II, researchers must have a “user friendly” method of obtaining and 
analyzing the necessary high resolution output data and validating the conjectures 
present in said algorithms. Additionally, as a starting point for the mathematical 


analysis, researchers must have a thorough understanding of the basic principles as- 


sociated with the various methods of calculating force attrition. This capability for 
analysis and foundation for understanding do not exist outside of the United States 
Army Concepts Analysis Agency (CAA) and its proprietary implementation of the 
COSAGE model and ATCAL Phase I and II algorithms. The intent of this thesis is 
to develop the framework for a test bed capable of conducting detailed analysis of the 
attrition problem. By developing realistic high resolution ground combat simulations 
using the Janus simulation software, capturing the detailed killer/victim results and 
the associated time series of casualties, and applying these results to the study of the 
aggregated attrition problem, continued and specialized future research can be con- 
ducted in this discipline. Specifically, as an immediate result of this research, future 
inquires can be made into the validity and appropriateness of the currently accepted 
versions of the ATCAL Phases I and II algorithms aiding in the proper documenta- 
tion of their conjectures and results. Moreover, with this vehicle for comparative 
analysis, efforts can expand to the development of alternative attrition algorithms to 


meet the needs of future aggregated modeling efforts. 


D. JANUS SIMULATION SOFTWARE 


As the initial step in creating a functional test bed for the study of attri- 
tion, one must choose a vehicle to conduct the sample data base of high resolution 
simulations. Janus is the primary high resolution, interactive wargaming simulation 
software used by the United States Army for the modeling of brigade sized (and be- 
low) operations, and, its widespread availability within the military community, ease 
of use, and verified credibility make it an obvious choice for meeting the demands of 
this proposed analysis. Named for the two-faced Roman god who was the guardian 
of portals and the patron of beginnings and endings, Janus had its beginnings as a 
nuclear effects modeling simulation developed by the Lawrence Livermore National 
Laboratory of the University of California. [Ref. 3] Its applicability to tactical train- 


ing and operations analysis quickly made it functional to the United States Army 


Training and Doctrine Command (TRADOC), Training Analysis Center (TRAC). 
The current version is intended to satisfy both the combat development and training 
communities. Its versatility allows it to be used interactively as a man-in-the-loop 
simulation to replicate the effects of command and control on a realistic battlefield 
during the simulation as enemy activity develops. As an analysis tool, Janus can 
be used non-interactively thus reducing the variability introduced by human decision 
making during the simulation execution. Janus is a stochastic simulation in that the 
system determines the results of actions (like direct fire engagements) according to 
laws of probability. While the principle focus of Janus is on ground maneuver and 
artillery units, it also models weather and its effects, day and night visibility, engineer 
support, minefield employment and breaching, rotary and fixed wing aircraft, resup- 
ply operations, and chemical environments. The inherent flexibility of this simulation 


software allows its user to accurately model a multitude of combat scenarios. 


E. SCOPE AND LIMITATIONS 


This research is unique in that there has been no other attempt to establish 
an independent capability for analyzing the aggregated attrition problem, specifically 
the ATCAL Phase I and II algorithms. This work will validate the capability of con- 
structing a high resolution ground combat simulation in Janus for the purpose of 
deriving attrition rate parameters using the ATCAL Phase I and II algorithms. Ad- 
ditionally, there has been no attempt to configure the Janus output data to conform 
to the requirements of the current ATCAL Phase I and II algorithms (currently exist- 
ing in the FORTRAN computer language). This procedure and the transformation of 
the output data is practical only within the confines of analysis and algorithm valida- 
tion because of the time expenditure required to develop the ground combat scenario 
in Janus and the nonexistence of automated software to link the Janus killer/victim 
results to current aggregated models. With initial efforts confined to the manual com- 


pilation of Janus output data and the explicit formulation and calculation of basic 
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mathematical models, this research could be expanded to focus on the compatibility 
and automation of the data linkage between Janus the current and/or future ATCAL 
computer codes. However, this is not the focus of this thesis. And, as only one specific 
ground combat simulation will be presented as part of this thesis, follow-on work is 
applicable in the development of additional scenarios with varying mission objectives, 
force definitions, environmental conditions, and terrain selections in order to general- 
ize results otherwise limited to specific simulation parameters. A robust collection of 
simulations would be required to properly validate the conjectures presented in the 
problem statement. As such, the results presented in this work relate only to the 
specific scenario and situations modeled and should not be considered to be the final 


solution. 


I] 





ae THE SIMULATION SCENARIO 


In order to generate the required high resolution output data with Janus for 
the development of attrition rate numbers using the ATCAL Phase I and II algo- 
rithms (analogous to the process described in Chapter I, Section B), the author has 
chosen to simulate a heavy brigade combat team (BCT) combat scenario. The typical 
heavy BCT includes a heterogeneous mix of direct fire, indirect fire, and supporting 
maneuver systems. Generically, the force would include tanks, armored personnel car- 
riers, artillery systems, air defense artillery systems, engineering equipment, wheeled 
vehicles, air assets, and obstacles. The BCT is the preferred level of combatant or- 
ganization for this study because of the described combined arms composition and 
the ability to portray this force in Janus. As a high resolution model, Janus can 
effectively support a force type of this size and capture the detailed results of all 
direct and indirect fire engagements. With these results, the attrition rate estimation 
process (ATCAL Phase I) and the deterministic battle replay (ATCAL Phase II) can 
be conducted. 

The tactical scenario chosen for this combat simulation is one of practical 
relevance. As this research focuses primarily on the modeling of ground combat 
as part of an overlying larger tactical operation, the author has chosen to simulate 
the mechanized task force operations in desert environmental conditions as a base- 
line scenario. As the Janus simulation package can be an interactive training tool 
when employed as a man-in-the-loop simulation, it is often used by TRADOC at the 
United States Army’s combat arms training centers. Specifically, the Armor Officer 
Advanced Course (AOAC) located at the Unites States Army Armor Center, Fort 
Knox, Kentucky, employs Janus as a tool to train and evaluate the mission planning 
and tactical decision making skills of its officer students. With a specified mission of 
training the current and future leaders of the United States maneuver warfare orga- 


nizations, AOAC and its focus on modern ground maneuver tactics, techniques, and 
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procedures provide an excellent simulation opportunity for this research. The Janus 
simulation used during this research was executed by AOAC Class 97-3 at the Fort 
Knox Brigade Battle Simulation Center. [Ref. 4] And, the simulation is available 
for review through coordination with Professor Bard K. Mansager, Department of 
Mathematics, Naval Postgraduate School. The stochastic nature of the simulation is 
defined not only by the probabilistic events which occur as a function of the intrinsic 
algorithms, but also by the human factors associated with the particular officer stu- 
dents who executed the simulation. If multiple runs of this scenario were executed, 


each would produce different results based upon these stochastic processes. 


A. MECHANICS OF THE SIMULATION 


At the highest level of organization, Janus allows the definition of six combat- 
ant sides where each side could represents an autonomous force. Under each side, 
forces can be further organized into task forces where modeled weapons and maneu- 
ver systems can be organized and manipulated in groupings that mirror their real life 
task organization. After the defined combatants are organized into task forces, one or 
more of these task forces are output to a specific computer workstation for user con- 
trol. Janus permits up to 24 computer workstations to be in use during a simulation. 
The 24 workstations are divided into side 1 through side 6 terminals based upon the 
number of users assigned to manipulate the respective sides and task forces. |Ref. 5] 
The simulation used for this research makes use of sides 1 and 2 and will hereafter 
be referred to as BLUE and RED, respectively. This assignment represents the roles 
of the combatants during the simulation. In this simulation, BLUE represents the 
friendly Unites States forces being controlled by the AOAC 97-3 student officers and 
RED represents a Soviet-style enemy forces being controlled by designated opposing 


force personnel. 
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1. BLUE Force Composition and Organization 

The BLUE force (BLUEFOR) is made up of the United States maneuver 
and weapons systems appearing in Appendix A, Table XIII on page 64; and, when 
organized for combat, these systems represents a heavy BCT. The BLUEFOR is 
considered to be “balanced” in terms of its armor to mechanized infantry ratio. As 
an integrated combined arms fighting force, the BLUEFOR systems are organized into 
a typical BCT hierarchical command and control structure as detailed by Table I. For 
the purposes of the AOAC student officers who executed the Janus simulation, this 
fictitious force was designated as the 1st Brigade Combat Team (1 BCT). 


Unit Description Major Systems 
4 Direct Fire Maneuver Battalions M1A1, M2 


2 Indirect Fire Artillery Battalions M109A6 
1 Air Defense Artillery Company BSFV 


1 Chemical Support Company IF TRU 
1 Engineer Support Battalion AVLM 

1 Logistics Support Battalion PRY 
1 Tactical Air Control Party A10 (6) 


Table I. BLUEFOR Task Organization for AOAC 97-3 Janus Simulation 
[Ref. 6] 


2. RED Force Composition and Organization 

The RED force (REDFOR) is made up of Soviet style equipment listed in 
Appendix A, Table XII on page 63; and, when organized for combat, it represents 
a less than full strength Soviet style motorized rifle regiment [MRR(-)]. Note the 
greater than three to one ratio of mechanized infantry to armor, implying that the 
REDFOR is better equipped for mechanized infantry operations than tank warfare. 
Additionally, the REDFOR has a significant number of indirect fire artillery systems 
at its disposal. The defensive REDFOR MRR(-) is organized within the hierarchical 
command and control structure detailed by Table II. For the purposes of the AOAC 
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student officers who executed the Janus simulation, this fictitious enemy force was 


designated as the 168th MRR(-). 


Unit Description Major Systems 


2 Direct Fire Maneuver Battalions BMP 

1 Direct Fire Tank Battalion T80 

4 Indirect Fire Artillery Battalions 251, 255, BM-21] 
1 Air Defense Artillery Company SA7, ZSU-23 


1 Engineer Support Company MTK-LI 
1 Logistics Support Company IF TRU 
1 Helicopter Company HIND (6) 


Table I]. REDFOR Task Organization for AOAC 97-3 Janus Simulation 
[Ref. 6] 


B. TERRAIN REPRESENTATION 

The Janus simulation uses digitized terrain developed by the United States 
National Imaging and Mapping Agency (NIMA). The terrain selected for this sim- 
ulation is located within the United States Army’s National Training Center (NTC) 
near Barstow, California, covering a 100 km? area with actual military grid coordi- 
nate locations defined by the corner points NK800400, NK790400, NK790390, and 
NK800390. The terrain is displayed in Janus similar to that used by military person- 
nel during actual training and combat operations. It includes contour lines (elevation 
information), roads, rivers, vegetation, and urban areas. The visibility issues as- 
sociated with terrain elevation and visual obstructions are realistically portrayed in 
Janus. These modeled attributes will degrade (or enhance) a combatant’s line-of-sight 
(LOS) and movement. The modeled vehicles and equipment must maneuver under 
the conditions represented by the digital terrain and its overlying natural and man- 
made features. For example, only those vehicles with a “swim” capability for crossing 
water obstacles are permitted to do so without employing bridging assets; and, when 


executed, the vehicles will proceed at a reduced movement rate. 
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The environmental conditions and time of battle are also controlled in Janus. 
The user can specify environmental conditions including temperature, cloud condi- 
tions, and air characteristics affecting smoke propagation. Weather and its effects 
influence the visibility and movement of the modeled combatants. The start time 
can be set in order to execute the combat simulation at a particular time of day. 
This accounts for simulating day and night operations which also affect visibility and 
movement attributes. The environmental conditions used during this simulation were 
summer, desert with 14 kilometers (km) visibility, and a 0800 start time. Since Janus 
accurately models the BLUEFOR and REDFOR weapons systems as a function of 
each system’s capabilities as defined in the database, users must consider all military 
aspects of employing their forces just as they would in actual combat. Invariably, 
Janus will reveal the weaknesses associated with poor planning and reward those who 
have tactically sound plans. As the combat scenario for this research was executed 
by students of the mounted warfare profession, the author assumes the employment 


of tactically sound techniques, procedures, and plans. 


C. SCENARIO EXECUTION SUMMARY 

The initial dispositions of BLUEFOR and REDFOR are depicted in Ap- 
pendix A, Figures 1 and 2 on pages 65 and 67, respectively. The BLUEFOR mission 
is to attack in zone to destroy or cause the withdrawal of the 168th MRR(-) secu- 
rity zone forces who are controlling the key terrain in and around grid coordinate 
NK620130 (Avawatz Mountain Passes) [Ref. 6]. The BLUEFOR executes a sup- 
porting attack in the north followed by a main attack in the south. The REDFOR 
conducts a defensive operation with an objective of stopping the penetration of the 
BLUEFOR and defending the significant terrain as mentioned above [Ref. 6]. The 
REDFOR executes a layered defense with two distinct defensive belts supported by 


deliberate minefields, obstacles, and barriers. 
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At the onset, the BLUEFOR is positioned to conduct offensive operations 
using two distinct avenues of approach (AA). For the purpose of this research, the 
northern and southern AAs will be named AA North and AA South, respectively. The 
main effort of the BLUEFOR attacks along AA South with the supporting attack 
along AA North. As a matter of mounted warfare tactics, the BLUEFOR leads 
with reconnaissance assets in order to locate and direct fires upon the REDFOR’s 
defensive assets. Following closely behind are the BLUEFOR’s initial maneuver and 
engineering assets. This battlefield organization will support direct fire engagements 
with the REDFOR upon their detection by the BLUEFOR reconnaissance elements. 
The BLUEFOR engineering elements are positioned forward to facilitate the clearing 
of obstacles and breaching of REDFOR minefields. Next in the BLUEFOR order of 
battle is the bulk of the fighting forces primarily consisting of tanks and armored 
personnel carriers. The indirect fire artillery assets follow in the BLUEFOR combat 
formation. They are positioned in an attempt to be within range of the REDFOR 
while maintaining reasonable security from REDFOR counter-artillery operations. 
Finally, the support oriented vehicles and equipment make up the rear of the attacking 
BLUEFOR. Positioned well to the rear and in support of the BLUEFOR, six A-10 
close air support (CAS) aircraft are available to be employed as a coordinated part 
of the attack. 

The REDFOR, while maintaining a layered defense supported by maneuver 
obstacles and minefields, also maintains a counterattack capability with reserve forces 
located to the rear of their main defensive area. These combatants can be employed 
to reinforce weakened defensive areas or to counter a penetrating attack made by the 
BLUEFOR. As detailed by Table II, the REDFOR is primarily a mechanized force 
consisting of armored personnel carriers. As such, the REDFOR’s tanks are employed 
judiciously and in positions where they can have the most significant effect on the at- 
tacking BLUEFOR. One of the most significant capabilities of the REDFOR is found 
in the number of supporting indirect fire artillery assets. The REDFOR maintains 
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both a close and long range indirect fire capability while keeping the majority of its 
artillery assets outside of the maximum range of the BLUEFOR’s artillery units. The 
REDFOR emphasizes the use of its indirect fire systems and their associated radar 
systems to locate and conduct counter-battery fire against the BLUEFOR. As the 
simulation progresses, the REDFOR attempts to halt the advance of the attacking 
BLUEFOR. 


1. Execution Summary: 0 to 30 Minutes 

The BLUEFOR begins offensive operations by employing its leading recon- 
naissance assets in AA North. Shortly thereafter, the BLUEFOR and REDFOR 
make contact with each other and exchange direct fires. At approximately 11 min- 
utes, the REDFOR long range artillery assets fire special munitions placing two area 
minefields each containing 480 anti-armor mines to thwart the advance of the attack- 
ing BLUEFOR in the north. As of yet, the BLUEFOR attacking along AA South 
have not made contact with the REDFOR defensive positions in that area. During 
and immediately following the BLUEFOR’s minefield breaching operations in the 
north, a heated direct fire battle takes place. Both sides employ their indirect fire 
artillery assets as their reconnaissance assets fix each other’s locations on the battle- 
feld. With 28 BLUEFOR and 9 REDFOR combatants killed during this 15 minute 
localized battle, the BLUEFOR attacking along AA North are able to overwhelm the 
REDFOR defensive forces and break through the first defensive line in the north. 
Shortly thereafter, at approximately 27 minutes, the REDFOR long range artillery 
assets reposition to alternate firing positions for security purposes. 

see Appendix A, Figure 3 on page 69 for a pictorial representation of the 
battlefield at the beginning of the scenario. Table III below details the force attrition 
at the completion of the 0 to 30 minute time period. See Appendix A, Figure 4 on 


page 71 for a pictorial representation of the battlefield at 30 minutes. 
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BLUEFOR Losses REDFOR Losses 


Direct Fire 34 10 
Indirect Fire 7 4 
Minefields 4 0 
Cumulative Losses at 30 Minutes | 45 14 


Table III. Force Losses for AOAC 97-3 Janus Simulation at 30 Minutes 


2. Execution Summary: 31 to 60 Minutes 

The BLUEFOR, who met with success in the north, proceed east-northeast 
on a high speed AA and are met with little resistance from the REDFOR. At ap- 
proximately 40 minutes, the attacking BLUEFOR moving along AA South come in 
contact with the first layer of REDFOR defenders in the south. The majority of 
the BLUEFOR following the two attacks begin movement north-northeast to take 
advantage of the successful penetration of the REDFOR’s northern defenses. The 
BLUEFOR in contact in the south encounter and begin breaching operations on 
three linear minefields positioned in front of the REDFOR. As BLUEFOR recon- 
naissance forces identify REDFOR targets to their front they direct indirect artillery 
fires against them resulting in only marginal success. At approximately 45 minutes, 
the remainder of the REDFOR’s defending forces in the first defensive layer engage 
the attacking BLUEFOR from AA South. The REDFOR defense in the south is 
much more formidable than that which was found in the north. The BLUEFOR has 
some difficulty identifying penetration points through the REDFOR’s defenses. The 
BLUEFOR suffers a large number of casualties during this battle. The main body of 
the BLUEFOR compresses into a more compact movement formation and increases 
its movement rate. 

Table IV details the force attrition at the completion of the 31 to 60 minute 
time period. See Appendix A, Figure 5 on page 73 for a pictorial representation of 


the battlefield at 60 minutes. 
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BLUEFOR Losses REDFOR Losses 


Direct Fire 51 14 
Indirect Fire 9 4 
Minefields 4 0 
Cumulative Losses at 60 Minutes 64 18 


Table IV. Force Losses for AOAC 97-3 Janus Simulation at 60 Minutes 


3. Execution Summary: 61 to 90 Minutes 

By 65 minutes, the BLUEFOR attacking in AA South are massed at the for- 
ward edge of the linear minefields and conducting passage through the cleared lanes. 
At the same time, six A-10 CAS aircraft enter an aerial holding area for possible 
employment against the REDFOR in the first defensive belt. The REDFOR has a 
small ambush or counter-attack force made up of two tanks and two armored person- 
nel carriers positioned at the extreme south-southeast of its first defensive layer in 
order to block the advancing BLUEFOR from using a high speed AA. However, the 
overwhelming BLUEFOR combat power moving through this area quickly engages 
and destroys this team. At this point, the BLUEFOR begins to reorient its forward 
movement to the north to take advantage of the previous success there. At approxi- 
mately 83 minutes, the BLUEFOR attacking to the north-northeast along AA North 
make contact with the second layer of REDFOR defenses. The BLUEFOR begins an 
indirect artillery barrage against these REDFOR. 

Table V details the force attrition at the completion of the 61 to 90 minute 
time period. See Appendix A, Figure 6 on page 75 for a pictorial representation of 
the battlefield at 90 minutes. 

BLUEFOR Losses REDFOR Losses 


Direct Fire 65 25 
Indirect Fire 10 5 
Minefields c 0 
Cumulative Losses at 90 Minutes 82 30 


Table V. Force Losses for AOAC 97-3 Janus Simulation at 90 Minutes 
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4. Execution Summary: 91 to 120 Minutes 

At approximately 93 minutes, the six A-10 CAS aircraft depart the aerial hold- 
ing area to attack the second layer of the REDFOR defenses from the southeast to 
the northwest. Shortly thereafter, at 100 minutes, two REDFOR radar controlled 
anti-aircraft systems engage and destroy three A-10 CAS aircraft. The remaining 
aircraft adjust their course and continue their attack on the defending REDFOR 
combatants. Prior to their egress, the remaining three A-10 CAS aircraft are de- 
stroyed by REDFOR shoulder-fired anti-aircraft missiles located in the northeast. By 
105 minutes, the main BLUEFOR attacking elements clear the linear minefields and 
being their movement north-northeast along the high speed AA. At 116 minutes, these 
attacking BLUEFOR make contact and exchange initial fires with the remainder of 
the REDFOR’s second layer defensive forces in the central and southern portion of 
the maneuver area. Heavy exchanges of artillery occur during the development of 
this high intensity battle. 

Table VI details the force attrition at the completion of the 91 to 120 minute 
time period. See Appendix A, Figure 7 on page 77 for a pictorial representation of 


the battlefield at 120 minutes. 


BLUEFOR Losses REDFOR Losses 


Direct Fire 92 35 
Indirect Fire 24 9 
Minefields We 0 
Cumulative Losses at 120 Minutes 128 44 


Table VI. Force Losses for AOAC 97-3 Janus Simulation at 120 Minutes 


5. Execution Summary: 121 to 150 Minutes 

During the battle within the second layer of the REDFOR defense, six ad- 
ditional A-10 CAS aircraft are called forward from their staging area to attack the 
defending forces. The majority of the six aircraft are highly successful until their 


destruction at 131 minutes by radar controlled anti-aircraft fires. Heavy exchanges of 
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direct fires continue to occur and the attacking BLUEFOR divides its formation into 
two distinct attacks, one to the north-northwest and another along the high speed 
AA to the east-northeast. 

Table VII details the force attrition at the completion of the 121 to 150 minute 
time period. See Appendix A, Figure 8 on page 79 for a pictorial representation of 


the battlefield at 150 minutes. 
BLUEFOR Losses REDFOR Losses 


Direct Fire 139 49 
Indirect Fire 32 9 
Minefields 22 0 
Cumulative Losses at 150 Minutes 193 58 


Table VII. Force Losses for AOAC 97-3 Janus Simulation at 150 Minutes 


6. Execution Summary: 151 to 180 Minutes 

The BLUEFOR attack to the north-northwest meets with a tank heavy REDFOR 
defensive element including some reserve forces coming from behind the second de- 
fensive layer. The BLUEFOR reacts by passing additional tank assets through to 
the front of its formation in an attempt to counter the REDFOR tank threat with 
like firepower. At 158 minutes, the BLUEFOR in the northwest directs indirect ar- 
tillery fires against the REDFOR with little effect. At 165 minutes, the main body 
of the BLUEFOR which was following the lead attacking forces moves farther for- 
ward to concentrate fires. The direct fire battle in the north-northwest becomes 
a battle of attrition with the BLUEFOR maintaining numerical superiority. How- 
ever, this translated into a target rich environment for the REDFOR. A successful 
penetration of the defense occurs and BLUEFOR continues their northwestern move- 
ment past the REDFOR’s second layer of defense. As the breakthrough occurs, the 
REDFOR in the south take advantage of the flanking shot opportunities against the 
passing BLUEFOR and are rewarded with a large number of kills. This time period 
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is marked by the BLUEFOR’s largest number of direct fire losses as compared to any 
other 30 minute time period. 

Table VIII details the force attrition at the completion of the 151 to 180 minute 
time period. See Appendix A, Figure 9 on page 81 for a pictorial representation of 
the battlefield at 180 minutes. 

BLUEFOR Losses REDFOR Losses 


Direct Fire 254 63 
Indirect Fire 44 9) 
Minefields Di 0 
Cumulative Losses at 180 Minutes 325) 72 


Table VIII. Force Losses for AOAC 97-3 Janus Simulation at 180 Minutes 


7. Execution Summary: 181 to 211 Minutes 

The BLUEFOR divides into two relatively equally sized formations with one 
attacking north-northeast and the second, south of the first, attacking along the high 
speed AA to the northeast. At 188 minutes, the REDFOR commits its reserves in 
the north to counter the BLUEFOR attack in that sector. After some success, the 
numerically superior BLUEFOR destroys the employed reserves. At 198 minutes, the 
BLUEFOR in the south completely clears the second defensive layer and all of its 
associated minefields. At 206 minutes, the REDFOR commits addition reserve forces 
against the BLUEFOR attacking in the south along the high speed AA in an attempt 
to thwart the advance. At 211 minutes, the combat simulation is terminated simulat- 
ing the achievement of one or more of the BLUEFOR commander’s ending conditions. 
Battle termination standards may include attrition to a specified force level, attrition 
to a specified number of remaining maneuver systems, and/or the accomplishment of 
stated objectives. This particular combat simulation was undoubtedly stopped due to 
a self-imposed time constraint as the simulation was executed by and for the benefit 


of AOAC Class 97-3. 
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Table IX details the force attrition at the completion of the simulation. See 
Appendix A, Figure 10 on page 83 for a pictorial representation of the battlefield at 


the completion of the simulation. 


BLUEFOR Losses REDFOR Losses 


Direct Fire 270 79 
Indirect Fire 45 9 
Minefields 29 0 
Cumulative Losses at Completion Minutes 344 88 


Table IX. Force Losses for AOAC 97-3 Janus Simulation at Completion 


D. UTILIZATION OF THE SIMULATION OUTPUT 


One of the most advantageous properties of the Janus simulation software is 
its ability to capture information about the execution of the combat scenario. The 
force loss totals detailed in the tables above are just samples of the information 
available to the Janus analyst. The Janus post execution output files which are 
important to the attrition rate analysis of this research include the direct fire report, 
coroner’s report, killer/victim scoreboard, and the artillery impacts report. While 
several other reports are available, these four reports will provide the time series of 
casualties data required for the development of attrition rate estimates applicable to 
the ATCAL Phase I algorithm. With the required simulation output data in hand, 
what remains to be developed are the mathematical constructs which serve as the 


basis for the ATCAL Phase I and II algorithms. 
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If. ATTRITION METHODOLOGIES FOR 
AGGREGATED MODELS 


A. BACKGROUND 

Prior to any detailed analysis of the attrition rate estimate process as it applies 
to the ATCAL Phase I and II algorithms (the primary goal of this research), it is 
important to establish a foundation with respect to the various methods for calculat- 
ing attrition and determining the outcomes of battles within aggregated simulations. 
In general, nearly all attrition computation methodologies can be categorized by two 
techniques: the Anderson (Force-Ratio [FR]) model and some type of Lanchester- 
type model. Integral to the second attrition approach are the Bonder-Ferrell method 
and the Combat Analysis (COMAN) maximum likelihood estimation (MLE) method 
for determining attrition coefficients for use in Lanchester-type equations (of which 
variants are used in the ATCAL algorithms). Although the COMAN MLE ap- 
proach is not specifically used within any current models, its applicability to the 
attrition modeling cannot be overstated. Each of these distinct attrition methodolo- 
gies maintains a relevance to modern aggregated modeling. For example (although 
not all-inclusive), the Anderson (FR) method is one of the two options for assessing 
attrition in TACWAR 4.1, the Bonder-Ferrell method is employed within the United 
States Army’s Vector In Commander (VIC) model, and the ATCAL methodology is 
used in CEM and ITEM. (TACWAR also has an option to use the ATCAL attrition 
methodology rather than FR.) One of the most interesting facts of this research is 
that ATCAL uses none of the attrition methodologies (explicitly) as they are de- 
scribed below. Rather, ATCAL uses an ad hoc method for attrition rate estimation. 
This is truly one of the mysteries of ATCAL for which this research attempts to 


establish the vehicle for further analysis. 
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B. THE ANDERSON (FR) BASED ATTRITION MODEL 

In this theory, a heterogeneous force (its overall strength) can be described 
by one scalar value representing the consolidated individual systems’ contributions 
to the force as a whole. Each system maintains an assigned value (score) relative 
to some common weighting criteria. Summing over all systems multiplied by their 
respective value (score) results in the generation of an overall firepower index (FPI) 
for the force being represented. As this FPI is a single scalar value, the resulting 
process is characterized as a homogeneous force ratio attrition model. The explicit 
representation of the heterogeneous force is replaced by this single scalar measure 
of combat power. This strategy can be expressed mathematically for a hypothetical 


force, X, as shown below in Equation III.1. 


Melle = Ss SG (Firepower Score) (III.1) 
ee — index each representing a unique weapon system type, total of m 
S* = value (score) of one weapon system of type : 
X; = number of weapon system of type z 


Supposing that the FPI of two opposing forces (X and Y) are represented by 


ae 
, gives a measure of 


Ocak 


relative combat power at the instant at which the values are calculated. The FRyy 


FPly and FPly, respectively, then the force ratio, FRxy = 


may be used to determine force missions, unit postures, casualty computations, and 
mission success. Specifically, the FRxy can be indexed against personnel loss curves 
to determine attrition to modeled forces. The loss curves employed in this process 
are generally derived from the results of historical battles. This fact is considered to 
be one the inherent strengths of the Anderson (FR) method. On the other hand, one 
of the perceived weaknesses of this method is associated with the determination of 


the individual weapon system values (scores) as discussed below. 
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1. Determination of Weapon System Value 

There is no general agreement in the modeling community on how to deter- 
mine the individual weapon system values (scores), S*. The fundamental methods 
for determining firepower score values include analysts’ perception of combat value 
(the naive estimate), historical combat performance, measures of weapon firepower 
and support platform characteristics, and, finally, measures of what a weapon can 
kill. This last method attempts to compute the value (score) for a weapon as being 


proportional to the rate at which it destroys the value of opposing enemy systems, 





and, as such, becomes a circularly defined problem. Moreover, this approach’s de- 
pendence on the “rate” at which weapons destroy opposing weapons in battle infers 
the reliance upon assumed or derived non-negative scalar attrition rates (one for each 
weapon system firing at every opposing weapon system for both forces) which account 
for the highly detailed engagement process described in Chapter I. The mathematics 
of this method, called the AntiPotential-Potential (APP) method, amounts to solv- 
ing for the individual weapon system scores through a linear system of equations (in 
matrix form) as detailed below by Equations III.2 and III.3. The solutions to these 
equations are found by solving an equivalent system of eigenvalue problems [Ref. 1], 
for which the development is beyond the scope of this discussion. However, as a prac- 
tical example, Appendix B details the application of the Anderson (FR) method and 
the subordinate APP process to hypothetical high resolution output data representa- 
tive of the AOAC 97-3 Janus simulation described in Chapter II. The calculations are 
sampled from the heterogeneous battle between selected BLUEFOR and REDFOR 
weapon systems. Note that this computational excursion does not account for every 
weapon system modeled in the Janus Simulation in an attempt to keep the associated 
eigenvalue problem numerically solvable in closed form without the aid of automated 


software. 
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Force X 
Force Y 


Cy 


pi SS, (111.2) 


j=l 


S (111.3) 
7—) 

er A ae uc | m distinct weapon systems 

hint Ween coe Yn} n distinct weapon systems 


value of one type z system in X force 

value of one type 7 system in Y force 

matrix of rates at which one X; system kills Y;systems 
matrix of rates at which one Y; system kills X;systems 
proportionality constant 


proportionality constant 


2. Advantages and Disadvantages 


The use of historical data, both in the generation of personnel and equipment 


loss curves and in the determination of weapon system values, generates a certain level 


of model believability and personal comfortableness with the Anderson (FR) method 


of attrition. Moreover, aggregated simulations using this method are credited with 


reproducing comparable results to the recorded outcomes of actual battles. This 


trait is often considered important when attempting to validate and verify a model’s 


usefulness. The desire to extend these result to the modeling of future battles is a 


task which generates a considerable amount of discussion. With respect to determin- 


ing weapon system values (scores), the APP method and its associated eigenvalue 


problem are an attractive procedure for calculating homogeneous force values and 


resulting force ratios for use in aggregated models because of their mathematical ease 


of computation (regardless of their appropriateness). 
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Conversely, there are a significant number of limitations inherent to this pro- 
cedure. A major criticism of this method is that the APP scoring method may be 
inconsistent with the personnel loss curves which are generally based upon historical 
data. Additionally, since Equation III.1 is additive across weapon systems, there is 
no ability to represent the interactions among complementary weapons on the battle- 
field. This is an especially serious shortcoming considering the nature of the current 
battlefield and the synergies of modern weapon systems. One must exercise caution 
when conducting comparative assessments of force mix or force balance since the 
scalars FPly and FPIy do not discriminate between balanced and unbalanced forces 
(in terms of size). An FPly = 100 and an FPIy = 100 are equal forces using this 
method regardless of the number of systems represented or their respective values. 
The reader will also note that Equation JII.1 is linear in the number of weapons X; 
of each weapon type z. Therefore, a model employing the Anderson (FR) method 
cannot be expected to effectively answer specific questions about the contributions of 
individual weapons. For example, the ramifications of an unmanageably large num- 
ber of a specific weapon may be intuitively obvious, but, the FR model will linearly 
increase the overall FPI without regard to practicality. It is the opinion of the author 
that the overwhelmingly mysterious characteristic of the Anderson (FR) method lies 
within the development of the individual weapon system values (scores). The reader 
may have surmised that in order to apply the Anderson (FR) method, the analyst is 
forced to arbitrarily assume values (scores) for each of the modeled weapon systems 
or to rely upon the killer/victim results of high resolution simulation to develop the 
attrition rates necessary to employ the APP method of solving for weapon system 
values (scores). Associated with using the high resolution simulation results to derive 
the kill rates within the APP method for increased believability, the reader should be 
aware of some of the inherent weaknesses which plague the APP eigenvalue problem. 
The score values derived from the eigenvalue problem are extremely dependent upon 


scenario characteristics and the number of weapons systems which translate into tar- 
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get engagement opportunities. Changes to the attrition rates (the key elements to 
the APP method) will typically cause all weapon system score values to change for 
both forces in ways that are hard to describe or justify. Consider another author’s 


evaluation of the APP method of determining weapon system values: 


The eigenvector score values sometimes change in ways that are hard 
to explain and that have been called paradoxical by some. For example, a 
shift in fire distribution designed to increase the kill rates of a higher value 
enemy target can sometimes reduce the total value index of the firing force. 
Whether this is, in fact, paradoxical depends on how deeply the relationships 
are followed. If there are few of the high value target, then shifting fire away 
from a lower value but more numerous enemy system might very well result 
in a lower total value being killed by the firing force. 

Other anomalies, however, are harder to explain. The numeric values 
of the scores are sometimes oversensitive to small changes in the input kill rate 
matrices... [Ref. 1] 

It is with this knowledge that the analyst is drawn to developing similar 


methodologies using nonlinear value equations, a formidable excursion which is be- 


yond the scope of this research. 


3. Typical Implementation of the Anderson (FR) 
Method 


A generic aggregated combat model, like TACWAR may employ the Ander- 
son (FR) method for attrition calculations in the following way. Assume that a hypo- 
thetical model represents opposing forces down to the division level of organization. 
The model would compute attrition for the divisions in a combat sector with an up- 
date cycle equal to the overall time cycle of the model, possible every 12 or 24 hours. 
Prior to the calculation of attrition, the simulation schedules several events to update 
the forces. These tasks may include movement of forces into and out of the combat 
sectors, allocation of reinforcements, determination of which force is the attacker for 
the cycle, determination of the defender’s posture for the cycle, and computation of 
the results of any scheduled air battles for the cycle. At a point during the execu- 


tion of a generic aggregated combat simulation where the conditions are met for two 
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opposing divisions to engage in battle, the Anderson (FR) method would be invoked 
to determine the respective forces’ FPIs. In a simplified sense, the resultant FRxy 
determines the outcome of the battle with respect to losses. After computing (mostly 
estimating) the parameters required for the Anderson (FR) and subordinate APP 
calculations, the resulting FRxy values are indexed against the appropriate person- 
nel casualty curve to determine overall attrition levels. (Unfortunately, indexing the 
FRxy against the casualty rate and equipment loss curves may not be consistent 
due to differing computational methods. The FRyy is derived using a linear model 
including the APP process, while the casualty rate and equipment loss curves may 
have been derived nonlinearly from historical data.) During each time cycle of the 
simulation (possibly 12 or 24 hours), the APP eigenvalue problem must be recursively 
solved to determine the updated weapon system scores and the resulting forces’ FPIs 
due to their dependence on the number of individual weapon systems remaining. This 
entire process 1s simply repeated as the simulation continues through additional time 


cycles. 


C. LANCHESTER-TYPE DIFFERENTIAL ATTRITION 
MODELS 


Prior to discussing the Bonder-Ferrell and COMAN MLE methods for esti- 
mating attrition rates, the reader must gain a basic understanding of the overlying 
framework to which the attrition rates are applied. Lanchester-type models origi- 
nated with the early works of F.W. Lanchester, a British engineer and inventor, in 
an attempt to quantify the advantages of concentrating fires in modern warfare. He 
hypothesized that ancient warfare was merely a collection of one-on-one duels where 
force size primarily determined battle outcomes. The resulting mathematical model 
representing this ancient form of warfare is referred to as Lanchester’s Linear Law 
or Firer-Target | Firer-Target (FT|FT) attrition. On the other hand, Lanchester 


supposed that modern warfare maintained the ability to concentrate the effects of 
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multiple weapons on surviving targets, thus creating a many-on-one condition. Ina 
hypothetical battle between forces X and Y, this condition would make the force X 
casualty rate proportional to the number of surviving firing systems in force Y. And, 
the resulting casualty-exchange ratio would be inversely proportional to the force 
ratio. This mathematical model for modern warfare conditions is termed Lanch- 
ester’s Square Law or Firer | Firer (F|F) attrition. These underlying Lanchester 
models, with various extensions and improvements, supply the basic framework to 


the ATCAL Phase I and II algorithms, which will be discussed in Chapter IV. 


1. lLanchester’s Linear Law (FT|FT Attrition) 


The primary assumption of ancient warfare is that any weapon system (soldier) 
not engaged in combat must wait until an enemy soldier becomes available before 
joining the battle. The assumption implies that the number of soldiers of force X put 
forth on the battlefield to engage force Y is determined by the number of soldiers of 
force Y present. Therefore, under the conditions of ancient warfare, there should be no 
advantage to concentrating forces. A particular soldier from a force is either engaged 
in a one-on-one duel, or not. Intuitively, each force size is a function of time, where 
the continuous real time variables z(t), y(t), and ¢ are approximations to their discrete 
analogs in an actual battle. The change of each force size with respect to time depends 
upon the number of firers engaged in combat and can be represented mathematically 
by the (homogeneous) Lanchester Linear Law (FT|FT), Equations III.4 and III.5, 


x d 
uivovaess aad 2 represent the change in the size of forces X and Y, respectively, 


dt dt 
over time. 
dx 
aos —az(t)y(t) with z(0) = zo (II1.4) 
d 
= —bz(t)y(t) with y(0) = yo (III.5) 
z(t) = size of force X at time t 
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size of force Y at time t 


y(t) 

a = attrition rate at which an individual Y firer 
kills X targets 

6 = attrition rate at which an individual X firer 


kills Y targets 


2. Lanchester’s Square Law (F|F Attrition) 

Extending the assumptions of ancient warfare to modern battles, one gains the 
attribute that many firers can engage a single target, thus achieving the advantages 
associated with concentration of fires. Two general cases of combat developed out 
of this study, namely, area fire and aimed fire. First, Lanchester assumed area fire 
to be a situation where a force spreads its fires over a general area occupied by the 
enemy rather that aiming directly at individual enemy weapon systems. Assuming 
that these area fires are uniformly distributed over the target area and that the 
area is independent of the number of targets, the Lanchester Linear Law (FT|FT) 
was assumed to hold true for area fires. Second, aimed fire assumes that individual 
targets are identified and attacked by one or more opposing firers. In this case, 
Lanchester assumed that the attrition of force X depended on how many force Y 
firers were directing aimed fires at force X. Mathematically, the basic homogeneous 
form of Lanchester’s Square Law (F|F) can be described by a system of differential 
equations, Equations III.6 and III.7, where end eae represent the change in the 


dt dt 


size of forces X and Y, respectively, over time. 


d 

= = —ay(t) with x(0) = xo (111.6) 

dy 

i; = —br(t) with y(0) = yo (UE?) 
z(t) = size of force X at time t 


y(t) = size of force Y at time t 
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attrition rate at which an individual Y firer 
kills X targets 
attrition rate at which an individual X firer 


kills Y targets 


The homogeneous representation of Lanchester’s Square Law (F|F) can be ex- 


tended to meet the needs of heterogeneous forces by accounting for all of the possible 


firer/target combinations in the simulated battle. The resulting heterogeneous for- 


mulation of Lanchester’s Square Law (F|F) is provided by the system of differential 


equations, Equations III.8 and IIJI.9, shown below. 


dz z 
dt 
dys 
dt 
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b;; 
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X force 
Y force 
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j=l 
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t=1 


size of force X, system 2, at time t 

size of force Y, system 7, at time t 

attrition rate at which Y; firers kill X; targets 
fraction of Y; firers allocated against X; targets 
attrition rate at which X; firers kill Y; targets 
fraction of X; firers allocated against Y; targets 
(re Coe 


{Y1, Y25 Y35--- Bua 


(III.8) 


(III.9) 


The reader should be aware that Equations IJI.8 and III.9 represent one pos- 


sible form of the Lanchester Square Law (F|F) equations. While the equations pre- 


sented represent work which was primarily accomplished during the 1960s, the overall 


methodology has evolved over time. In the 1970s, fire allocation factors which influ- 


36 


ence the attrition rates were developed by submodels executed in support of the 
overall process. Nevertheless, the equations described above appear in a fundamental 


form which aids in the understanding of this topic. 


3. The Significance of a;; and };; 

As Lanchester’s original purpose was a qualitative justification of the principle 
of concentration as an advantageous tactic in modern warfare, the derivation of the 
heterogeneous attrition rate coefficients a;; and b;; was of no particular concern. How- 
ever, today, as Lanchester based attrition models are employed as part of theater level 
aggregated models, the significance of these attrition rates cannot be overemphasized. 
As part of the model verification process, current aggregated models must provide 
reasonable and quantitative results. These numerical values assigned to a;; and 6,; 
are what represent the synergies of battle between various weapons systems and the 
very nature of specific actions on the battlefield. Their derivation is what brings 
believability to the equations present within combat models. Several methods have 
been proposed and utilized over the years to derive the attrition coefficients. The 
remaining background discussion of this chapter focuses on two of these methods of 


attrition rate estimation which are currently in use. 


D. THE BONDER-FERRELL METHOD FOR ATTRITION 
RATE ESTIMATION 


The Bonder-Ferrell method for estimating attrition rates makes use of an in- 
dependent analytical model which does not depend on input from other sources. 
Instead, a stochastic process model is developed for each Y; firing at a single passive 
type X; target, and the expected time between casualties (E[T;*]) is determined, where 
T;, is a random variable representing the time it takes a single Y; firing system to 


kill one X; target. Assuming that the casualty process is a Markov (Poisson) process 


nd 


(see [Ref. 7]), the attrition rate parameters can then be expressed as a;; = 


eee 
E\T;; 
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j 
pute E[T] by calculating the distribution of the number of shots required to achieve 


y= . While Bonder’s original work focused on deriving the equation to com- 
a specified number of hits, Ferrell’s extension of this work modeled this process as a 
Markov renewal process (see [Ref. 7]) with three discrete states (new engagement, hit, 
and miss) and the associated transition probabilities and times. Thus, the expected 
time required to kill an enemy target is the long run expected time required within 
the Markov renewal process until moving on to a new engagement. As early studies 
in this area focussed on homogeneous forces, continued work in the 1970s extended 
the concept to heterogeneous battles by combining conditional kill rates with target 
priority lists. Additionally, target acquisition was characterised as either parallel or 
serial. 

One of the most significant advantages to the Bonder-Ferrell method of at- 
trition coefficient estimation is that the computations are based upon what is going 
on in the battle at the given moment. There is no requirement to reference a data 
base library of high resolution simulation output to match up the conditions of battle 
between the overlying model and the high resolution battle. Additionally, since the 
underlying data used in these computations are basically measurable, the results are 
more credible to the user. This method is generally considered to be restrictive in 
that the in-depth engagement model is analytic, thus suppressing some detail. This 
form of attrition rate estimation also precludes the ability to capture synergistic ef- 
fects of weapon system combinations as the underlying process does not account for 
weapon system interactions. The reliance upon engineering level data does lend this 
approach to the explicit definitions of engagements; however, the amount of data 
required on each weapon system is great. Relying upon closed form equations for the 
development of the a;; and 6;;, the Bonder-Ferrell method allows for the re-evaluation 
of the attrition coefficients at each time step of the simulation by simply recalling and 


recomputing the required input values. 
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EK. THE COMAN MAXIMUM LIKELIHOOD ESTIMA- 
TION (MLE) APPROACH TO ATTRITION RATE 
ESTIMATION 


The COMAN MLE attrition estimation methodology is a fitted parameter 
model which examines the time series of casualty times (from a high resolution sim- 
ulation) and computes the maximum likelihood estimate for the mean time between 
casualties. While other statistical methods exist for determining point estimates of 
this type (including the method of moments, Bayesian estimation, and the method of 
least squares), the MLE approach is the only method which has had any significant 
application in combat simulations [Ref. 8]. Similar to the Bonder-Ferrell approach, 
the COMAN MLE method relies upon the assumption that the time between casu- 
alties is a Markov (Poisson) process with the memoryless property (time until the 
next casualty does not depend upon the previous time before the last casualty). The 
COMAN MLE method is most easily explained by detailing its mathematical con- 
struct. The underlying assumption of this method is that the high resolution simula- 
tion data represents a sample from the continuous analogue to the deterministic homo- 
geneous equations of Lanchester’s Square Law (F|F) (see Equations [II.6 and III.7), 
1.e., the continuous-parameter Markov chain with transition probabilities given by 


Equations III.10 and III.11 [Ref. 9]. 


Prob[X casualty in a small time interval of length At] = an(t)At 
(III.10) 


Prob[Y casualty in a small time interval of length At] = bm(t)At 
(III.11) 
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Additionally, consider the following random variable definitions and their as- 
sociated properties required for the development of the MLE attrition coefficient es- 
timation process. The reader can assume that any realization of a described random 


variables will be notationally represented by its corresponding lower-case letter. 


1 ifthe k” casualty is an X combatant, 


ox = 
Q otherwise, 
ay 1 ifthe k** casualty is an Y combatant, 
k — 
Q otherwise, 
K 
G = Se c* total number of X casualties 
k=] 
K 
cr = > > c total number of Y casualties 
pa 


ee + ch = K total number of X and Y casualties 


n; the realization of the number of Y combatants 


pms) 
o~ 
o~ 
a 
Nee” 

I 


just after the occurrence of the k*” casualty 


(t,) my, the realization of the number of X combatants 
Mlk 


just after the occurrence of the k** casualty 


Interpreting the random variable definitions given above, note that during 
the time interval [t,,t,4;] there are m, X combatants and n, Y combatants alive 
for values of k = (0,1,2,...,A —1). Armed with high resolution simulation output 
data of the form {(ti,c*,cY ), (t2,c*,ct),...,(tx,c%,cxr)}, the statistical estimates 
for the continuous time Markov chain representation of the deterministic form of 
Lanchester’s Square Law (F|F) can be computed using the method of maximum 
likelihood estimation (MLE). As a battle develops, the casualty processes for the X 
and Y forces are simply two superimposed Poisson processes with respect to time (a 
random variable). Developing and maximizing these likelihood functions from their 


applicable joint probability density function (PDF) (representing the times between 
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casualties) will result in the derivation of the attrition rate estimates, @ and 6. The 
maximum likelihood estimation process can be summarized by the steps detailed in 


Table X. 


Step 1 | Find the probability density function for the time until an X casualty. 
Find the probability density function for the time until a Y casualty. 
Construct the respective likelihood functions. 


Step 2 
on“ (i.e., the density functions for the observed sequences of events) 
Step 3 | Determine the values of a and b that maximize the 
respective likelihood functions. | 


Table X. Maximum Likelihood Estimation Process Steps 

















[Ref. 9] 
1. Maximum Likelihood Estimation Steps 
Applying the three maximum likelihood estimation steps to the stochastic rep- 
resentation of Lanchester’s Square Law (F|F) represented by Equations IJJ.10 and HI.11 
will result in the statistical derivation of 4 and 6. The author has chosen to detail 
the development of the estimated attrition parameters @ and b from the X casualty 
point of view. Initially, C and 1 — C~ will be used to represent X and Y casualty 
occurrences, respectively. 
a. Step 1 
Let the random variable Sy denote the time between X casualties with 
the well-known exponential PDF fs,(s-) = anexp ‘*")**, Similarly, the PDF for 
the time between Y casualties is fs,(sy) = bmexp~°™)*¥. Further, let the random 
variable S denote the time between any two consecutive casualties be defined as the 
minimum of Sx and Sy. The resulting PDF describing the time between any two 


casualties is 


fs(s) = (an + bm) exp7@7t8™s | (11.12) 


4] 


Incorporating the previously defined variable Cj‘, which equals one when the k** 
casualty is an X combatant and zero when the k” casualty is a Y combatant, one 


can further specify Equation IJJ.12 above. Note that Cf’ = 1 if Sy < Sy. The 


probability of this occurrence is Pwolell( OR =) = Seis Likewise, a Y casualty 
an + bm , 
occurs if Sy > Sy which implies Cf = 0 and Prob(C~ = 0) = a Therefore, 
an + bm 


the joint probability of the random variables S and Ci is 


Pro <3. Gas = 0) (1 —exp7©"t!™)*) . Prob(c* = 2) 


where 0 << s < coand:=0,1 


an x bm anes : (1 _ go ee 


(— Se | 


an + bm 


an+bm 


Therefore, the representative joint PDF of S and C7 is 


an x bm aoe Mins 
Ts,ox (s, cp ) = grin k. agg oe oe (an at bm) exp (an+bm) 
= (an) %* (bm)'~% exp-(ertim)s (111.13) 


where c* = 0,1. 


b. Step 2 

The likelihood function is developed by considering that the K casu- 
alties have occurred at times (t),t2,t3,...,tx) with a total number of X and Y ca- 
sualties represented by c{ and ch, respectively, where c* + ch = K. Recalling that 
C= 1 — 6, — 0 and C2 — 0 —— ©, — Bauation fin trcan nesmnctiona iy, 


written as 


(arg_y)* (brrg_1)% exp (re-1FomMa=1 (tet) | (III.14) 
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Equation III.14 represents the contribution of the k** casualty to the likelihood func- 
tion. As the times between casualties are independent random variables, the likeli- 
hood function for the observed sequence of events, Equation III.15, is the product of 


all of the independent contributions. 


K 
L(a,b) = [] (are—i)* (bray_1)* exp  (2re-1 thm a1 (te te-1) 
k=1 
(11.15) 
Gi Step 3 


Finally, the COMAN MLE process requires the maximization of this 
likelihood function. However, it is a common mathematical practice to maximize the 
natural logarithm (In) of the likelihood function as detailed below by Equation III.16. 
Both maximum values will occur at the same point [Ref. 10] and the logarithmic form 


is more easily manipulated. 


K K K 
In[Z(a,6)]) = © cg In(ang_1) + D> cj In(bme_1) — D> (ang—a + bmu-1) (te — te-1) 
k=1 k=1 k=! 
(II1.16) 
Continuing with Step 3, one can take the partial derivatives of Equation III.16 with 


respect to a and 6 to get Equations III.17 and III.18, respectively. 


Oln{L(a, b)] eS 
a ae do ne—-i(te — te-1) (11.17) 
y kK 
i nh Ce (III.18) 
Ob b k=1 


Finally, maximize Equations IJI.17 and III.18 by setting each of the partial deriva- 
tives equal to zero and solving for a and b. The resulting solutions for a and 8, 
denoted @ and 6, maximize the likelihood function (Equation III.15) and serve as the 


COMAN MLE estimates for the attrition rate coefficients. See Equations IIJ.19 and IJI.20. 


op 
i ee eee III.19 
es ye CT iets aD, 

. Y 
a (III.20) 


yo mp-i(ty — te-1) 


2. Summary of the COMAN MLE Method 

While the presented mathematical description of the COMAN MLE method 
for estimating attrition rate coefficients is applicable only to the homogeneous rep- 
resentation of Lanchester’s Square Law (F|F) in one of its simpler forms, the reader 
can see that the general method is quite simplistic. One of the simplifying charac- 
teristics of the developed example is the appearance of the attrition rate coefficients 
as linear terms only. It is quite possible that the true nature of the attrition rate 
coefficients is non-linear in actual combat; however, for the purposes of this dis- 
cussion, the author has chosen to represent the mathematics of the method with 
tractable forms. As a practical example, Appendix C details the application of the 
COMAN MLE method for attrition rate estimation to the high resolution output data 
from the AOAC 97-3 Janus Simulation described in Chapter II. The calculations are 
presented for the mechanized battle between the BLUEFOR M2 (with TOW missile) 
and REDFOR BMP-2 (with AT-5 missile). Note that this computational excursion 


makes use of Lanchester’s Square Law (F|F) (as the mechanized battle is one char- 


44 


acterized by exchanges of aimed fire) and the mathematical framework leading up to 


and including Equations III.19 and III.20. 


F. EXTENSION OF ATTRITION RATE CONCEPTS 
TO THE ATCAL ALGORITHMS 


By no means should the reader consider this treatment of attrition rate anal- 
ysis to be exhaustive. A significant number of noteworthy modifications have been 
proposed and are used within Lanchester-type deterministic models. For a highly 
detailed inspection of these modifications, the interested reader should consult the 
prolific works of Professor James G. Taylor, Naval Postgraduate School, and co- 
advisor to this research. However, with the general background development of the 
methods for attrition rate estimation complete, the reader is now prepared to consider 
their application to (or lack thereof) the ATCAL algorithms. The specific modifica- 
tions to the underlying Lanchester-type model which form the basis for the ATCAL 
methodology can be primarily attributed to Gordon M. Clark and his PhD Thesis, 
“The Combat Analysis Model” |Ref. 8], while ongoing analysis is being conducted by 
James G. Taylor [Ref. 2]. However, recall that the specific attrition rate parameter 
estimation process used in ATCAL Phase I can be described as an ad hoc process at 


best. 
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sh THE UNDERLYING 
LANCHESTER-TYPE MODEL FOR ATCAL 


As the primary purpose of this research is to establish a test mechanism for 
evaluating the current ATCAL attrition algorithm as it applies to aggregated combat 
models (like CEM and TACWAR), the interested reader must develop a familiarity 
with ATCAL’s underlying Lanchester-type models and their application via explicit 
assessment equations. As detailed documentation describing the theoretical bases of 
ATCAL is not available from its original author, this chapter attempts to outline the 
construct of the ATCAL attrition algorithm based primarily upon a very thorough 
pe ent of the subject by James G. Taylor. [Ref. 2] It is with an understanding of 
the current ATCAL algorithm that the interested reader could conduct research for 
improving its mathematical efficiency or developing mathematical variants to fulfill 
other combat model requirements. The reader may find it helpful to review Chap- 
ter III, Section C, which details the basic framework of Lanchester’s Linear (FT|FT) 
and Square (F|F) Laws because of the parallelisms which extend to the Lanchester- 


type models present in the ATCAL attrition algorithm. 


A. THE HOMOGENEOUS FORCE MODEL 

Considering the case of homogeneous combat between two opposing forces, the 
ATCAL attrition algorithm employs the nonlinear Equations IV.1 and IV.2 which can 
be attributed to Gordon M. Clark [Ref. 8]. 


= = —a(1—p™)y(t) with #(0) = zo el) 

d 

= = —B(1—q%)z(t) with (0) = yo (IV.2) 
HG) 0 size of force X at timet 
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yi) 0 size of force Y at timet 


ee 0 rate at which an individual Y firer kills acquired X targets 
Geo rate at which an individual X firer kills acquired Y targets 
ee 1 probability that an individual Y firer has no 


X target available for engagement 
0< gl <] probability that an individual X firer has no 


Y target available for engagement 


The reader may find it more convenient to express Equations IV.1 and IV.2 
with respect to target availabilities rather than non availabilities. ‘To do so, let 1—(1— 
A)*() = 1 — p*), which is the probability that an individual Y firer has one or more 
X target(s) available for engagement. And, analogously, let 1—(1— B)# = 1—q¥. 
With these substitutions, Equations [V.1 and IV.2 can be expressed as follows. 


dz 


= = ~afl— (1 - A) ly (2) with 2(0) = xo (IV.3) 
4 = —6[1 — (1 — B)* z(t) with y(0) = yo (IV.4) 


Considering the elements of Equation IV.3, the reader will note that the mul- 
tiplication of the conditional kill rate (a) with the probability of target availability 
(1 — (1 — A)*")) produces the expected rate at which an individual Y firer kills X 
targets. Similarly, the multiplication of @ and 1—(1— B)# yields the expected rate 
at which an individual X firer kills Y targets. 

It has been shown that force levels can never become zero or negative as a 
matter of theory (see [Ref. 9]). However, as the attrition calculations are necessarily 
executed on computer platforms, the potential for achieving “computer zero” must 
be recognized. 5o, in practice, the attrition models above must be defined to be valid 


only when the forces levels, z(t) and y(t) are greater than zero. This refinement of the 
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homogeneous attrition model within the ATCAL algorithm produces Equations IV.5 
and IV.6. 


Sill = AAU IGTEA ES tore eel vc) 
de _ f -all—(1—A)*]y(t) for x(2), y(t) — 
dt 0 otherwise. 

—B{1—(1 — B)\¥) z(t) for z(t), y(t) > 0, 
dy _ [ BU - (= Bye) for a(t), v(e = 
dt | 0 otherwise. 


This positive force level requirement can be easily incorporated in computer algo- 
rithms by including a conditional check of the current force level minus any calcu- 
lated change to the force level. If the resulting force level is less than zero, then the 
algorithm should return a value of zero. Otherwise, the current force level is updated 


by subtracting the calculated attrition. 


1. Application of the Averaging Operator 
In order to transform Equations I1V.5 and IV.6 into explicit assessment equa- 
tions as employed in the ATCAL algorithm, it is important to develop the mathe- 


matical properties of the averaging operator (over time) given by 


Ave|z|—= 3 = -[ r(s)ds (IV.7) 


where z(t) is a function of time. 
Application of this operator to the Lanchester-type equations developed above will 
result in an operator equation for the force levels, z(t) and y(t), respectively. Define 
the number of casualties inflicted on the X force as Az and the number of casualties 


inflicted on the Y force as Ay. 
Ar = zo—z(t), where zp = x(0) 
Ay = yo— y(t), where yo = y(0) 


, aoe , 
Now, consider applying the averaging operator to Ti which results in 


dx l srtdz 
Avel—| = -—-|/ —- 
oer t Jo ds : 
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= sIz(t) — Zo] 

_ _[to-2(%)| 
t 

7" Ax 

a 


Hence, applying the averaging operator to the homogeneous Lanchester-type 
Equations IV.5 and IV.6 representative of the ATCAL algorithm, one obtains the 


operator equations, Equations IV.8 and IV.9 below. 


S = —ofl —(1— A) ly(2) 
= Ave = Ave[—all — (1 — A)" ]y(2)] 
ee SF = —aAvelll — (1 - 4) le) 
<=> Ar = atAvel[1— (1 — A)]y(t)] (IV.8) 
and, similarly, 
Ay = ftAvel[1— (1 — B)Y]z(t)). (IV.9) 


The reader will note that the arguments of the averaging function appearing in 
Equations IV.8 and IV.9 above are highly nonlinear. Unfortunately, these equations 
cannot be further simplified to explicit expressions involving only average force levels 
without making some type of approximation. The key approximation made within the 
ATCAL algorithm is that the average of this nonlinear argument in z(t) and y(t) is 
equal to the same functional relationship between z and y, the results of the averaging 
function applied to z(t) and y(t). Although this approximation is not precise, it is 
often applied to simplify averages of nonlinear functional relationships. It is this 
simplification that allows ATCAL to move forward computationally. Equations IV.8 
and IV.9 can now be replaced by Equations IV.10 and IV.11. 


Az = at[{1—(1—A)*]y (IV.10) 
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Ay = @t[l1—(1- B)*)]z (IV.11) 


Zz = _ time average of the force level z(t) 

y = time average of the force level y(t) 
Ax = number of X casualties during time t 
Ay = _ number of Y casualties during time t 


2. Applying the Results for Homogeneous Forces 

As written, Equations [V.10 and IV.11 have four unknown quantities, namely, 
the average force levels (% and y) and the numbers of casualties (Az and Ay). With 
two equations and four unknowns, additional assumptions are required in order to 
solve this system of assessment equations. The second key assumption made in the 


ATCAL algorithm is that average force levels can be represented by the expressions 

(—Az) __ _{—Ay) 
jens OOO rag 
in(1 — 27) ade) 
hold true, respectively, if z(t) and y(t) are exponential functions (increasing or de- 


i . [Ref. 11] It so happens that these expressions 
creasing) over time. As force levels necessarily decrease during battle, the negative 
exponential function is assumed. Combining this information with Equations IV.10 
and IV.11, the underlying Lanchester-type attrition model used to determine the 
numbers of casualties during each assessment period is represented in the ATCAL 
algorithm by the following system of four simultaneous nonlinear equations in the 


four unknowns, z,y, Az, and Ay. 


RI 

| 

—] | 
| |p 
aID8 
W 


nad 

lI 
= 
7 
ye 


While small perturbations of this notation are explicitly modeled in the ATCAL 
algorithm, this system of equations should be thought of as the basic ATCAL casualty 


ol 


assessment system of equations. This form better suits future research in applying 
numerical solution techniques to solve the system (recognizing the requirement for 
excluding negative force levels as described in Equations IV.5 and IV.6). The actual 
method used in ATCAL to solve this system of nonlinear equations is a heuristic 
form of successive substitutions which is beyond the scope of this research. The 
interested reader may find the study of alternative numerical solution techniques 
extremely relevant to verifying current ATCAL methodologies and proposing future 


replacement ones. 


B. THE HETEROGENEOUS FORCE MODEL 

Considering the case of heterogeneous combat between two opposing forces 
made up of m distinct X force weapon systems and n distinct Y force weapon systems, 
the ATCAL attrition algorithm employs the nonlinear Clark Equations IV.12 and 
IV.13 detailed below. In these equations, 2 = 1 denotes the lowest priority target 
type and 2 = m denotes the highest priority target type within the X force. And, 
7 = 1 denotes the lowest priority target type and 7 = n denotes the highest priority 
target type within the Y force. Again, the reader may find it helpful to compare 
these equations with Equations III.8 and III.9 of Chapter III. 





da; “ : Z | 
i = —>>a,;(1 - py”) II pty, Font — 2) eae 
= k=i+] 
(IV.12) 
dy, bi 
a = — >) By — oF oS) I qe 23 oye Fl, Doe yn 
= my 
(IV.13) 
et) 0 number of weapon systems of type z in force X at time t 
DAT number of weapon systems of type 7 in force Y at timet 


(0) = Pepe tie eal rae 


o2 


y;(0) = Y;0 105 Jom 2,-4.7 


Cie eal rate at which an individual Y; firer kills acquired X; targets 
0 rate at which an individual X; firer kills acquired Y; targets 
0< pi <1 probability that an individual Y; firer has no 


X; target available for engagement 


Ws qj ) =< hh probability that an individual X; firer has no 


2 


Y; target available for engagement 


Similar to the homogeneous case, the reader may find it more convenient to express 
Equations IV.12 and IV.13 in terms of target availabilities by letting 1 -—(1— Ajj) = 
1- pa” and 1 — (1— B,,;)¥") =1- ja With this substitution, Equations IV.12 


and IV.13 can be expressed as follows. 


—— co ENC. = Ae) sil eet ee) t= 1,2,---,m 
dt j=l k=i+1 

(TV .14) 
dy; = =e a _ (1 = B;,)¥@™) II (1 -_ By)* }2;(t) We lane AU 
dt =I k=j+1 

(IV.15) 


Again, in order to preclude the possibility of achieving “machine zero” (or 
negative force levels) during numerical attrition calculations, Equations IV.14 and 
IV.15 can be further specified so that the force levels z;(t) and y;(t) are valid only 
when they are greater than zero. One final notational simplification renders the 
heterogeneous formulation of the underlying Lanchester-type model for ATCAL more 
compact. Consider the set of indices corresponding to target types with a higher 


priority than a given target type by formally defining 


I;; = {kl|k € integers} 


X target type k has a higher priority than target type z for a Y; firer; and, 
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Jy ab exe tniegers) 


Y target type & has a higher priority than target type 7 for a X; firer. 


With this simplifying set notation, a preferred form for the underlying Lanchester- 


type equations for ATCAL appear below in Equations IV.16 and IV.17. This rela- 


tively simple representation results in less computational requirements when imple- 


mented as the attrition model of a large-scale aggregated combat model. 


dix — DL a {(1 — we) Maen, Pej }ys(t) for 2i(t),ys(4) > 0, 


dt 0 


0 


otherwise. 


(IV.16) 


otherwise. 


djy — 17, Bf — a) Mees, FE Fait) for 2i(t), y(t) > 0, 


(IV.17) 


1. Application of the Averaging Operator 


Without loss of generality, the application of the averaging operator (with 


associated notational definitions) to the heterogeneous Lanchester-type equations will 


result in operator equations for the force levels, z;(¢) and y;(t). In order to maintain 


consistent notation, the author has chosen to represent the equations using the “target 


availability” notation of Equations IV.14 and IV.15. 








a aisf( —(1— Ajj)7) Ul (1 — Axs)** }y;(t) 


k=1+1 


Avel-Y aiif( — (1 — Aj;)**) ll (1 — Agj)*) Fy; (t)] 


k=1+1 
— VawAvel{(1— (1 Ag) TT 0 — Au) ns (0 
j=1 =i+1 
>» aaj Ave[{(1 — (1 — Aig)**) TI (1 — Ags)** }y;(t)] 
j=1 =i+1 


(IV.18) 


o4 


and, similarly, 


Tr 


Ay; = £3- BjrAve[{(1 — (1 — B,)”) Nee ee) ja): 


t=] k=9+1 


(IV.19) 


Analogous to the homogeneous case, the key approximation made within the 
ATCAL algorithm is that the average of these nonlinear arguments in z;(¢) and y;(t) 
are equal to the same functional relationships between Zz; and y;, which are the results 
of the averaging function applied to z;(t) and y;(t). Applying this approximation, 
Equations [V.18 and IV.19 can be replaced by Equations IV.20 and IV.21. 


Az = tS ais {(1 — (1 - Ay)*) I (1 — Ans)"* $9; (IV.20) 
3=1 k=i+1 


td BiK(1 = ea) ll a ere. (ZA) 


k=j+1 


Ay; 


2. Applying the Results for Heterogeneous Forces 

Similar to the homogeneous argument presented in the previous section, both 
the average force levels (Z; and y;)and the numbers of casualties (Az; and Ay;) 
are unknowns in Equations IV.20 and IV.21. Therefore, this system is made up of 
(m+n) nonlinear equations in 2-(m-+n) unknowns. In order to increase the number 


of equations to equal the number of unknowns , ATCAL again makes use of an 


wa! . —Az; 
assumption for the average force levels resulting in the expressions 7; = a 
‘ ee =) 
and y; = ASAY i [Ref. 11] Combining this information with Equations IV.20 


In(1 — a2) 
and IV.21, the underlying Lanchester-type attrition model used to determine the 
numbers of casualties during each assessment period is represented in the ATCAL 
algorithm as the following system of 2- (m+ 7) simultaneous nonlinear equations in 


the 2- (m+n) unknowns, Z;,y;, Az;, and Ay;. 
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Ati get Wiar t = (1 = Aaj) 7) a eso 
Ay; = 1072; Sat —  — 8;,)”) Wi eee ee, 
a m8 ) 
7. eeu) 
2 amas) 


While small perturbations of this notation are explicitly modeled in the ATCAL 
algorithm, this system of equations should be thought of as the basic ATCAL caualty 
assessment equations. This form better suits future research in applying numerical 
solution techniques to solve the system (recognizing the requirement for excluding 
negative force levels as described in Equations IV.16 and IV.17). The actual method 
used in ATCAL to solve this system of nonlinear equations is a heuristic form of 
successive substitutions which is beyond the scope of this research. The interested 
reader may find the study of alternative numerical solution techniques extremely 
relevant to verifying current ATCAL methodologies and proposing future replacement 


ones. 


C. SUMMARY OF THE LANCHESTER-TYPE ASSESS- 
MENT EQUATIONS IN ATCAL 


The homogeneous and heterogeneous systems of nonlinear equations found 
at the ends of Sections A and B, respectively, comprehensively describe the under- 
lying Lanchester-type models and their basic assessment equations for the ATCAL 
attrition methodology. “The most significant implication of these results, however, 
is that a much more convenient way of generating numerical results for ATCAL 
has immediately been established: namely, numerical integration of the underlying 
Lanchester-type differential equations.” [Ref. 2] Moreover, with these mathematical 
constructs, the interested reader is now prepared to further investigate the soundness 
of the heuristic successive substitution method of ATCAL and compare it to more 


robust numerical solution techniques. 
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V. CONCLUSION 


A. CONCLUSIONS 

This research has developed a practical test mechanism for studying the attri- 
tion methodologies of aggregated combat models. The Janus high resolution simula- 
tion software, capable of modeling brigade sized (and below) operations, is completely 
adequate for producing the required high resolution output data required for the anal- 
ysis of aggregated attrition methodologies. Janus is both widely available within the 
military community (to include the Naval Postgraduate School) and is easily learned 
through the use of a tutorial. See [Ref. 5]. This user friendly computer simulation 
Biavides the researcher with sufficient detail and a mechanism for capturing detailed 
interactions on the modeled battlefield. The ability to capture the highly detailed 
interactions between weapon systems is what qualifies Janus as a valuable tool for 
conducting analysis of the attrition algorithms inherent to ATCAL (and other ag- 
gregated attrition methodologies). The reader should follow the steps outlined in 


Section C to initiate such a study. 


B. RECOMMENDATIONS 

First, having established the vehicle (Janus) for generating the data required 
for developing attrition (kill) rate parameters, the author has opened the door for 
continued research in the area of attrition methodologies. As this research developed 
and discussed one such combat scenario, the continued study of attrition calculations 
would require additional Janus scenarios with varying parameters. The creation and 
execution of such scenarios could easily be incorporated into current high resolution 
modeling courses taught at the Naval Postgraduate School. 

second, as discussed in Chapter IV, the nonlinear attrition assessment equa- 
tions resulting from the Lanchester-type equations used in the ATCAL algorithm are 


currently “solved” using a heuristic method of successive substitutions. Significant 


of 


follow-on research is applicable in the employment of a more robust numerical so- 
lution technique for these systems of equations. This area of continued study does 
not necessarily depend upon additional Janus scenarios; however, additional scenarios 
would provide a data source for technique verification. 

Third, the computer code which contains the ATCAL Phase I and II algo- 
rithms exists in several mediums. A significant endeavor integral to making ATCAL 
(and any proposed modifications) more efficient includes a thorough analysis of the 
current computer code which exists in at least two languages, FORTRAN and C. 
The interested reader may find that the algorithms, which make up ATCAL would 
be more efficient in one or more specific computer languages. For example, a proposed 
numerical solution technique may be more conducive to one computer language than 
another. 

Finally, an automated bridge is required between the Janus high resolution 
output data and the computer code used to perform the aggregated attrition calcula- 
tions (like ATCAL). While this capability may not be directly applicable to the actual 
execution of aggregated models as they currently exist, it is integral to the efficient 
study of the algorithms themselves. It is in the research arena that this automated 
capability for linking Janus output to attrition algorithms would be helpful. JETS 
may be an important tool for continued study of this recommendation. 

As combat models are an integral part of the military analysis and training 
communities, continued analysis of attrition methodologies (specifically ATCAL) in 
the field of aggregated combat modeling will greatly enhance model believability and 
overall usefulness. As computer technologies continue to outpace attrition algorithm 
development, verification, and validation, a continued focus towards efficiency and 


accuracy will ensure the applicability of current and future attrition methodologies. 
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C. INITIATION OF ATTRITION ANALYSIS USING 
JANUS 


The key to initiating a detailed study of attrition and its inherent algorithms 
lies within the generation of high resolution combat simulation data. As such, the 
researcher must carefully consider the process of creating high resolution scenarios. 
The interested researcher should consider the following steps which outline the model 


formulation and data collection process. 





Define the tactical scenario (forces composition, terrain type, 
environmental conditions, and concept of the operation) 
for which a model is desired. 


Step 2 | Construct the Janus simulation scenario and validate its 
pi? | pont ofthe deine aeration between deed forte 
Step 3 | Determine an appropriate number (7) of iterations for which to 
| |e Jann cimsaton, en 
Step 4 | Independently execute the Janus simulation n times to develop 
a robust sample of simulation output data from which statistically 
significant killer/victim information can be determined. 
Step 5 Analyze the Pre high resolution output data and estimate 
rr i ation Ql rater bet opening tepateeyarts 
Step 6 | Incorporate - ao attrition (kill) rates in the analysis of 


attrition algorithms inherent to pat (and other aggregated 
attrition methodologies). 


Step 7 | Utilize an appropriate Lanchester-type deterministic model to 
replay the simulated battle. 


Table XI. Steps for Attrition Analysis and Deterministic Battle Replay Using Janus 














59 





rr 


LIST OF REFERENCES 


[1] James K. Hartman, Sam H. Parry, and William J. Caldwell. Aggregated Combat 
Modeling. Naval Postgraduate School, Monterey, California, 1992. 


[2] James G. Taylor. Working Papers, SAC of OSD PA & E Project. 1997-1998. 


[3] Titan. Applications Group Inc. The Janus 3.X/UNIX Model User’s Manual. 
Fort Leavenworth, Kansas, 1993. 


[4] Robert Gordon. Janus Brigade Simulation, AOAC 97-3 (software). Fort Knox, 
Kentucky, 1997. 


[5] Bard K. Mansager, Gerald M. Pearman, Larry R. Larimer, and Michael L. Shenk. 
Janus Version 6.0 Tutorial. Number NPS-MA-97-007. December 5 1997. 


[6] Small Group 2 Armor Officer Advanced Course. 52d Infantry Division (Mecha- 
nized) Operations Order 97-06-01. June 1 1997. 


[7] Sheldon M. Ross. Introduction to Probability Models. Academic Press, Inc., San 
Diego, California, 1993. 


[8] Gordon M. Clark. The Combat Analysis Model. Columbus, Ohio, 1969. 
[9] James G. Taylor. Lanchester Models of Warfare, Volume IJ. March 1983. 


[10] Richard J. Larsen and Morris L. Marx. An Introduction to Mathematical Statis- 
tics and Its Applications. Prentice-Hall, Englewood Cliffs, New Jersey, 1986. 


[11] U.S. Army Concepts Analysis Agency (CAA). ATCAL: An Attrition Model 
Using Calibrated Parameters. Technical Paper 83-3, Bethesda, Maryland, 1983. 


[12] Gilbert Strang. Linear Algebra and Its Applications. Harcourt Brace & Company, 
Orlando, Florida, 1988. 


[13] William E. Boyce and Richard C. DiPrima. Elementary Differential Equations 
and Boundary Value Problems. John Wiley & Sons, Inc., New York, New York, 
1992. 


61 





APPENDIX A. JANUS COMBATANTS 
LISTING AND RUN-TIME BATTLEFIELD 
PICTURES 


This Appendix provides the reader with the detailed description of the combat 
systems making up the REDFOR and BLUEFOR modeled in the AOAC 97-3 Janus 
Simulation described in Chapter II. Additionally, pictures are included to depict the 
battlefield details and force allocations prior to the start of hostilities and at 30 minute 


increments throughout the simulation. 


T-80- Main Battle Tank 15 
BMP-2 Armored Personnel Carrier 55 
BRDM-A Armored Personnel Carrier 10 
BRDM-2 Armored Personnel Carrier 7 
GMZ-MI Armored Mine Warfare System 2 


MTK_LI Armored Engineer System ] 
2512 M 122mm Self-Propelled Artillery 16 
IF TRU Wheeled Truck 26 
DF Tru Wheeled Truck 16 
251 155mm Self-Propelled Artillery 30 
BM-21 Rocket Launcher 8 
SNAR-2 Radar 2 


ZOOPAR Radar 2 
255 152mm Self-Propelled Artillery 3 
ADA TM Air Defense Artillery Team 5 
ZSU-23 Air Defense Artillery 5 
HIND Attack Helicopter 6 
POL AI Aviation Fuel Depot 2 


Total REDFOR Combatants 242 
Table XII. REDFOR Organization for AOAC 97-3 Janus Simulation 
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M1A1 Main Battle Tank 82 


M1A1P Main Battle Tank w/ Mine Plow 24 
M1A1R Main Battle Tank w/ Mine Roller 8 
M2 Armed Armored Personnel Carrier ae 
BSFV Armed Improved Armored Personnel Carrier 6 
M113 Armored Personnel Carrier 21 
FIST-V Artillery Command Vehicle 4 
CEV Combat Engineer Vehicle 4 
BRIDGE Mobile Bridge Carrier 5 
ADA TE Air Defense Artillery Team 4 
ACE Armored Combat Earthmover 13 
AVLM Armored Vehicle Launched Bridge 9 
AVLM-W Armored Vehicle Launched Bridge 6 
AVLB W Armored Vehicle Launched Bridge 2 
HMV.50 Wheeled Vehicle 47 
HMV/M1 Wheeled Vehicle 10 
HMV/SM Wheeled Vehicle 6 
HMV/CO Wheeled Vehicle 8 
IF TRU Wheeled Truck Das, 
DF Tru Wheeled Truck 4 
M106A1 120mm Mortar Vehicles 22 
M577 Armored Personnel Carrier (Command) 6 
GSR/M1 Ground Surveillance Radar 1 
SCT TE Scout Team 34 
SAPPER Combat Engineer Team 90 
SAW TE Machine Gun Team 144 
VOLCANO Mine Dispensing System 4 
M109A3 155mm Self-Propelled Artillery 24 
M109A6 155mm Self-Propelled Artillery 18 
HMMQ36 Artillery Radar System 1 
ADA TE Air Defense Artillery Team 4 
SMK113 Armored Vehicle Smoke Generator 3 
A-10 Fixed Wing Close Air Support 12 
Total BLUEFOR Combatants 764 


Table XII]. BLUEFOR Organization for AOAC 97-3 Janus Simulation 
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Figure 4. Battlefield Picture at 30 Minutes for AOAC 97-3 Janus Simulation 
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Figure 5. Battlefield Picture at 60 Minutes for AOAC 97-3 Janus Simulation 
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APPENDIX B. APPLICATION OF THE 
ANDERSON FR METHOD TO A 
HYPOTHETICAL HETEROGENEOUS 
BATTLE REPRESENTATIVE OF THE 
AOAC 97-3 JANUS SIMULATION 


This appendix provides the reader with an example application of the Ander- 
son FR Method with subordinate APP process for attrition rate estimation using 
hypothetical kill rate data representative of the simulation output of the AOAC 97-3 
Janus Simulation described in Chapter II. Specifically, the author has selectively 
chosen to calculate the BLUEFOR and REDFOR Firepower Scores (FPS) and the 
resultant Firepower Index (FPI) using two unique BLUEFOR weapon systems and 
three unique REDFOR weapon systems. The heterogeneous aimed fire battle between 
these opposing weapons systems provides the reader with a numerically obvious closed 


form representation for both forces’ FPS. 


1. FORCES COMPOSITION AND NOTATIONAL TRANS- 
LATION 


Let X represent the BLUEFOR composed of two weapon system types, namely, 
114 M1 Main Battle Tanks and 113 M2 Bradley Fighting Vehicles. Let Y represent 
the REDFOR composed of three weapon system types, namely, 15 T-80 Main Battle 
Tanks, 55 BMP-2 Armored Personnel Carriers, and 17 BRDM Armored Personnel 
Carriers. Recall from Chapter III that the goal of the Anderson FR Method is the 
representation of the heterogeneous force as a single scalar measure of combat power. 
Translating Equation III.1 from Chapter III, we have the mathematical expressions 


detailed below by Equations B.1 and B.2. 
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2 
bel = pes (BLUEFOR Firepower Score) (B.1) 


ol 


1 = integer index where 1 =M1 and 2 =M2 


S* = value (score) of one M1 
Ss = value (score) of one M2 
X1 eel | evils 


X42 = [3 Nes 


3 
make = Se (REDFOR Firepower Score) (B.2) 
j=1 


j = integer index where 1 =T-80, 2 =BMP-2, and 3 =BRDM 
SY = value (score) of one T-80 
SY = value (score) of one BMP-2 
SY = value (score) of one BRDM 
Yast; 
Y2 = 55 BMP-2s 
¥3 = 17 BRDMs 


2. DETERMINATION OF WEAPON SYSTEM VAL- 
UES 
Recall that the fundamental methods for determining the values of S* and 
on include the analysts’ perception of combat value, historical combat performance, 
measures of weapon firepower and support platform characteristics, and, measures of 


what a weapon can kill. The last method attempts to compute the value (score) for 


a weapon as being proportional to the rate at which it destroys the value of opposing 





enemy systems. This method for determining weapons system scores, called the 


AntiPotential-Potential (APP) method, is the primary focus of this appendix. 
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The formulation of the APP problem relies upon the “rate” at which systems 
destroy opposing systems in battle. These assumed or derived non-negative scalar at- 
trition rates (one for each weapon system firing at every opposing weapon system for 
both forces) account for the highly detailed engagement process described in Chap- 
ter I. While the actual kill rates for each pairwise combination of weapon systems 
resulting from the Janus simulation described in Chapter II could be derived through 
analyzing the killer/victim scoreboard, the author has chosen to use hypothetical kill 
rate values for each pair of interacting weapon systems thus defining kill rate ma- 
trices, K;; and Lj; below. This simplification assists the reader with focussing on 
the task at hand, understanding the Anderson FR methodology and the subordinate 
APP process. An example derivation of attrition rates from high resolution output 


is discussed in detail in Appendix C. 


a 
i. Ae 
Jase = and Le: ee 
eo 
2 8.3 

The matrix K;; represents the rate at which one type 2 weapon system (1=M1 and 
2=M2) kills type 7 weapon systems (1=T-80, 2=BMP-2, and 3=BRDM). Analo- 
gously, L;; represents the rate at which one type 7 weapon system kills type 2 weapon 
systems. As an illustration, the reader will note that the BLUEFOR (X force) M1 
Main Battle Tank kills the REDFOR (Y force) BMP-2 Armored Personnel Carrier 
at the rate of K1. = .5. 

The remaining unknowns of Equations B.1 and B.2 are the weapon systems 


scores, S* and SF. Deriving these individual scores equates to solving Equations B.3 


and B.4 (below) for the unknowns C'x, Cy, S;*, and SY. 


3 

OS SS fori = 1,2 (B.3) 
oe 
2 

Cae — SLs for 7 = 1,2,3 (B.4) 


t=1 
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Cx 


proportionality constant 


Cy = proportionality constant 


So, for example, the value of the BLUEFOR M1 Main Battle Tank (S*) is propor- 
tional to (C'x) the sum total value of every REDFOR system that it kills. The reader 
will note that the first equation above represents two equations (one each for 7 = 1, 2) 
and the second equation represents three equations (one each for 7 = 1,2,3) for a 
total of five equations. Moreover, this system of five equations has seven unknown 
variables, Cy, Cy, S*, S#,S*,SY, and SY. A system of linear equations with more 


unknowns than equations defining them is referred to as underdetermined. 
a. Developing the Representative Eigenvalue Prob- 
lem 


To solve the value equations defined by Equations B.3 and B.4, consider the 


following substitutions. 


2 
aS rota eel 


i 


aS lee Sepa ror 4) deed 


a 


Cy Si 


a 


Substituting the expression for Sy above into Equation B.3 yields 


ns = a LiS*)] tore) SN 


3 

d Ki 

ii 21 
2 3 
=> CxCyS§ = DISK LiS") for: =1,2, (B.5) 


and, applying the same construct to Equation B.4 


oe 2 
Cx Cys — Sapa o for 7 = 1,2,3. (B.6) 


j=12=]1 
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For clarity, Equations B.5 and B.6 can be expressed more compactly us- 
ing vector and matrix notation as depicted by Equations B.7 and B.8. Note that 
2 1 (Kj Li) = (KL) is a square matrix of dimension (2x2) and 4_, 1 (LyiKij) = 


4 ee re | 
(LK) is a square matrix of dimension (3x3). 


2600 3500 .4000 
4100 4 


2500 .3400 


ece= and LK =| .2000 .2700 .3100 


.1400 .1900 .2200 





GC se 
COM 


(KL)S* (B.7) 


(LK)S* (B.8) 


By letting AX = C'xCy [Ref. 9], the reader will note that Equations B.7 and B.8 
can be written as depicted below and represent a pair of eigenvalue problems which 


will share common eigenvalues. 


NSS = WGLyS* (B.9) 
\SY = (LK)SY (B.10) 


The reader should note that the matrix products (KL) and (LK) are square 
matrices which have only positive entries. Consequently, the Perron-Frobenius The- 
orem states that the largest eigenvalue of a matrix with only positive entries and 


the entries of the eigenvector associated with the largest eigenvalue are positive real 


numbers. See [Ref. 12]. 


b. Solving the Representative Eigenvalue Problem 
The eigenvalue problems represented by Equations B.9 and B.10 can be solved 


by finding the roots of the respective characteristic equations, 


det(KL—2I) = 0 


89 


4100—X  .5500 
<=> det( = 0 
2500 .3400 —» 
(Sani 
and, 
det(LK—XI) = 0 
2600 —»  .3500 4000 
<> det( 2000 00) =e O00 ee 0 
1400 1900 .2200—A 


where det is defined as the determinant and / is the identity matrix of the appropriate 
dimension. These characteristic equations are simplified to the following equations 
by evaluating the determinant and collecting like terms. If required, the reader can 
review the methods for evaluating determinants of (2x2) and (3x3) matrices in [Ref. 


12). 


det(KL—XI) = (.4100— 2)(.3400 — A) — (.2500)(.5500) 
= ? —.7500A + .0019 
and 
det(LK — I) = (.2600 — A){(.2700 — )(.2200 — 2) — (.1900)(.3100)} 


—.3500[(.2000)(.2200 — A) — (.1400)(.3100)| 
+.4000{(.2000)(.1900) — (.1400)(.2700 — )] 
=> 4.750002 00a 


= —X(d? —.7500A + .0019) 


Setting both determinants equal to zero and solving for the values of A which 


fulfill these equalities (using the quadratic formula [Ref. 12]), we find that the first 
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characteristic equation yields roots Ay = .7475 and Az = .0025, and, the second 
equation yields A, = .7475, A2 = .0025, and A3 = 0. The reader is now able to confirm 
the existence of the largest eigenvalue, A* = .7475, which is both real and positive 
as ensured by the Perron-Frobenius Theroem. Substituting this value, A* = .7475, 
into Equations B.9 and B.10 presents the following eigenvector problems for which 


we seek the solutions. 


| 
Cl 


CRI = 0 (B.12) 


l 
C1 


Cre 1S” (B.13) 


Consider Equation B.12 in isolation. 


_ 


PL 1S a0 
—.3375  .5500 ce 0 
—— 
.2500 —.4075 Se 0 


As there are infinitely many eigenvectors which solve this equality, the author will 


arbitrarily fix S* = 1 to produce a unique value for S¥* (and hence, a unique eigen- 
vector). This amounts to normalizing the value of the BLUEFOR M1 Main Battle 
Tank equal to one, which will force all of the weapon systems values to be scaled 
relative to the M1. The reader should be aware that this scaling method is only one 
from a list of many methods. Further information can be found in [Ref. 9]. Solving 


for the second unknown (the value of the BLUEFOR M2 Bradley Fighting Vehicle), 
l 


6134 | 


However, the reader should be aware of the fact that this solution vector is 


the reader can confirm that SX = .6134. Therefore, S¥ = 


not unique. Any scalar multiple of this value vector will also satisfy Equation B.12. 
The proportionality constants, Cx and Cy, appearing in Equations B.7 and B.8 will 
simply adjust to absorb the scalar multiple changes to the eigenvector of weapon 


system values - a condition for which the chosen scaling method shortly corrects. 
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Now consider the effectiveness of the BLUEFOR M1 Main Battle Tanks as 

they destroy the “value” of the REDFOR weapon systems. This can be evaluated 
3 

mathematically by the expression See os . The IDAGAM scaling method [Ref. 


j=l 


3 
9] seeks to scale the values of Sr such that 5 > hans = V\*, where \* = .7475 
= 


(the largest eigenvalue). Recall that in the development of the eigenvalue problems 

represented by Equations B.9 and B.10 that the author let 41 = CyCy. Combining 
3 

Equation B.3 with the known information S$ = 1 and >> ky;S; — V/)* the reader 


j=1 
will notice the following result. 


3 


ORS OSG SE 


=1 
Cx = VF 
Cx = /CxCy 
(Cx)? = CxC 
Cy = Cy =vV-* 


3 
Substituting the kill rate information from matrix K into the equation )> K 1553 = 
j=1 


VA* results in the equation .4S7 + .59Y + .5S¥ = .8646. This particular scaling 
method now calls for the replacement the first equation of the linear system repre- 
sented by Equation B.13 with this new equation. This action results in the following 


matrix equation. 


(Hig= 7 1)\S6 = 20 
4000 .5000 5000 on .8646 
<=> | 2000 —.4775" 31100 Sy 0 
1400 .1900 —.5275 oS 0 
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Several methods exist for solving this system of three linear equations in three un- 


knowns (many of which are detailed in [Ref. 12] and [Ref. 13]). When solved, the 
8176 


resulting values for the REDFOR weapon systems are SY = 6309 
4442 


3. CALCULATING THE RESULTANT FIREPOWER 
INDICES 


Armed with the values for each of the BLUEFOR and REDFOR distinct 
weapon systems, the reader can easily calculate the single-valued scalar descriptors 


FPIy and FPly. 


2 
Relea Xs: 
oS 
= (1)(114) + (.6134)(113) 
= 183:31142 
3 
Bivesey Sg 
j=1 


= (.8176)(15) + (.6309)(55) + (.4442)(17) 
= 54.5151 


As described in Chapter III, the Force Ratio (FRxy) can be calculated as FRxy = 
FPIx — 183.3142 


FPIly 54.5151 


termine the outcome of a battle with respect to parameters set within an aggregated 





= 3.3626. In a simplified sense, this resultant FRxy would de- 


model. For a more detailed treatment of the theoretical foundations of the Anderson 
(FR) Method and a critique of its use in force-ratio-based attrition models, see [Ref. 
9]. 

The reader may have gathered that the complexity associated with hetero- 


geneous battles made up of many distinct weapon systems would require the intro- 
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duction of an automated capability to solve the representative eigenvalue problems. 
For example, thoroughly accounting for the 33 BLUEFOR and 18 REDFOR distinct 
weapon systems modeled in the Janus simulation described in Chapter II would result 
in a linear system of 51 equations in 53 unknowns (counting Cx and Cy). The reader 
may also surmise that the kill rate matrices associated with this complex simulation 
may be sparse; in other words, many of the entries of the kill rate matrices could be 
equivalent to zero based upon weapon systems which do not directly interact on the 
battlefield. This fact introduces another consideration to the Anderson FR method 
and the subordinate APP process — the possibility for ill-conditioned matrices (see 
[Ref. 12]). At this point, the reader should have a basic understanding of the An- 
derson (FR) based attrition model and its application to a simplified heterogeneous 


battle. 
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APPENDIX C. APPLICATION OF THE 
COMAN MAXIMUM LIKELIHOOD 
ESTIMATION METHOD TO A 
HOMOGENEOUS BATTLE WITHIN THE 
AOAC 97-3 JANUS SIMULATION 


This appendix provides the reader with an example application of the COMAN 
MLE method for attrition rate estimation using portions of the simulation output 
data from the AOAC 97-3 Janus Simulation described in Chapter II. Specifically, the 
author has chosen to apply this method to the time series of attrition between two 
opposing mechanized infantry weapon systems. The aimed fire battle between the 
BLUEFOR M2 Bradley Fighting Vehicles and the REDFOR BMP-2 Armored Per- 
sonnel Carriers is examined as a homogeneous exchange of direct fires. The primary 
weapon systems engaged in this battle are the TOW anti-tank missile (BLUEFOR) 
and the AT-5 anti-tank missile (REDFOR). As such the derived attrition rate pa- 
rameters will be applicable to the Lanchester Square Law (F|F) for aimed fire war- 
fare. Finally, the author develops a practical application of the differential equations 
(Equations C.1 and C.2 detailed below) to validate the derived attrition rate esti- 
mates. This replay is analogous to the replay executed within the ATCAL Phase II 
algorithm. (However, the reader should recall that ATCAL Phase I uses an ad hoc 
method for attrition rate estimation rather than the COMAN MLE method.) The 
data for this analysis was exported from the Janus simulation output file and manipu- 
lated in a commercial spreadsheet software package (not included in this Appendix). 
The data appears in Table XIV on page 104 of this Appendix. However, the in- 
terested reader can obtain copies of the associated spreadsheet files from Professor 


Bard K. Mansager, Department of Mathematics, Naval Postgraduate School. 


95 


ile NOTATION 

Let X represent the BLUEFOR M2 systems (with TOW missiles) and Y 
the REDFOR BMP-2 systems (with AT-5 missiles). The relative Lanchester Square 
Law equations, modified with the COMAN MLE attrition rate estimates, model the 


change in force size with respect to time. These are Equations C.1 and C.2. 


d 
= = —ay(t) with z(0) = zo (C.1) 
dy ° 
os —bz(t) with y(0) = yo (C.2) 
z(t) = size of force X at timet 
y(t) = size of force Y at timet 
a = derived attrition rate at which BMP-2 firers kill M2 targets 


= derived attrition rate at which M2 firers kill BMP-2 targets 


Recall Equations III.21 and III.22 on page 44 of Chapter III, as shown below. 


oi 


re mee tra) 


A> 


er 


ee Tp oe mae bana) 


2, DEVELOPMENT OF THE ATTRITION RATE ESTI- 
MATES 


Applying Equations [II.21 and III.22 to the Janus output data found in Ta- 
ble XIV, we have the following: 
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ee 
per M1 (te — te-1) 
Total No. of M2 Casualties (Kk) 
y,(No. of BMP-2 Surviving at t, — 1) - (Time between k:, and k — 1,, Casualty) 
17 
469493 seconds 
0.00217256 M2 kills per BMP-2 firer per minute; and, 


A 
Sa myp—1(tp — te-1) 
Total No. of BMP-2 Casualties (K’) 
yo (No. of M2 Surviving at t, — 1)- (Time between k,, and k — 1,, Casualty) 
12 
997883 seconds 
0.00072153 BMP-2 kills per M2 firer per minute. 


Intuitively, from the standpoint of classic military operations, the resulting 


attrition rate estimates are logical with respect to the modeled combat scenario. The 


REDFOR is executing a layered defense with fortified positions and established fields 


of fire, while the BLUEFOR is on the attack characterized by weapon system exposure 


and less distinct fields of fire. Therefore, as one would expect, the REDFOR BMP-2 
kills the BLUEFOR M2 at a higher rate than vice versa (all else being equal). This 


provides the analyst with a common sense oriented way of justifying his or her inter- 


mediate results. 


DERIVATION OF FORCE LEVEL EQUATIONS FROM 
LANCHESTER’S SQUARE LAW (F|F) EQUATIONS 


From Equations C.1 and C.2, the informed reader may see that these ordinary 


differential equations can be solved to obtain expressions for force levels as a function 


of time. Consider the following differential equations: 
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dz 


Pig —ay(t) 
dy ~ 
a = —bz(t) 


which, when differentiated with respect to the independent time variable t, yield 


ch , dy 
de ~ "dt 
TA a 
TSE dt 


ee d d ne vw 
Substituting the values for - and 7 and using the initial conditions asso- 
ciated with Equations C.1 and C.2, one has the following two (second order, linear, 
constant coefficient, and homogeneous) ordinary differential equations (ODE). For a 


more detailed explanation of ODEs of this type, consult [Ref. 13]. 





dr a3 
TOE —ir5 initial condition 1 
Ez) = —dYo initial condition 2 
dt 
d*y * 
qe = bay(t) (C.4) 
y(0) = yo initial condition 1 
a —bro initial condition 2 


As such, each second order ODE can be solved using a characteristic equation 
after assuming general exponential solutions of the form z(t) = c; exp” +c. exp™’ 


r3t 


and y(t) = c3 exp™’ +c, exp”. Considering the second order ODE Equation C.3 and 


its initial conditions in isolation, suppose that z(t) = exp”, where r is a parameter 
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dx si —— 
to be determined. Then it follows that or rexp”™ and aT r? exp”. Substituting 


these values into Equation C.3 yields 


r*exp’’ = dabexp”™ 


(r? —ab)exp’ = 0. 


And, since exp” > 0 for all values of t, then (r? — 4b) must equal zero. Therefore, 
r = +1/ab are the roots to the characteristic equation. Note that @ and b are both 
positive real numbers implying that Tee are both positive real numbers. Knowing 
this, the complex and repeated root cases can be eliminated [Ref. 13]. General 


solutions to the second order ODE (Equation C.3) can now be represented by 


20) = cy exp’ #8)! 4 65 exp (Vat, 


Examining the initial condition 1, 


ci) cg 
<7 C1 +C2 = Zo, 
and initial condition 2, 
dzx(0) 





dt 
<3 Veh -aVab = —aYyo. 


The reader will recognize that the resulting two equations representing the 
initial conditions applied to our general solution form a system of two linear equa- 


tions in two unknowns, namely, c; and cz. Solving this system of equations yields 


l a l a pene: 
c= 5 (=. _ i] and c) = 5 (= +. (i). Substituting these values back into 


Og 


the general form of the solutions, Equations C.5 and C.6 result. (Without loss of gen- 
erality, Equation C.6 follows from the work presented for Equation C.5; the interested 


reader is encouraged to verify the result.) 


x(t) ; (=. — so] exp(V ab - (= - ve) op Vn (C.5) 
; ( = /$x) exp(V ar ap ( a ($2) op Van (C.6) 
a a 


As an alternative to the exponential representation of the force size equations 
exp’ + exp? 
2 


y(t) 


(dependent on timet), one can make use of the hyperbolic functions cosh = 


exp” =x exp? 


and sinh = to transform Equations C.5 and C.6 into 


8 
—— 

~, 
—— 

| 


ro osh | (V/ab) 1 -w 5 sinh {(Va8) f (C2) 
ip 708 (vai) | a aE sital (vai) J | (C.8) 


The most significant advantage to this form is the insight provided for the parametric 


.S 

oo 
o~ 

—’ 
| 


relationships that exist within the force level equations themselves. Consider the 
comparison of force size at time t compared to the initial force size for a hypothetical 


force X. This can be accomplished by dividing Equation C.7 by zo, resulting in 


=~ (8) ]-Byp lO) eo 


From Equation C.9, the reader can see that the level of force X appears to be de- 
pendent upon three terms; the initial force ratio (*) , the relative fire effectiveness 
ro 


measure (=) , and what can be called the intensity of combat indicator ( ab). In 


a battle to annihilation, these indicator values can be used to generalize simulation 


results [Ref. 9]. 
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Hence, one can now deterministically solve for force sizes, x(t) and y(t), for 


the desired time ¢ using Equations C.5 and C.6 (or Equations C.7 and C.8) and the 


estimated attrition rate coefficients @ and b. This practical extension of Lanchester’s 


Square Law (F|F) affords the analyst a straightforward method for estimating force 


size throughout a battle. The force size results can also be used to verify a and b 


themselves by comparing the numerical values of x(t) and y(t) to the analogous 


high resolution results (force size) at various points throughout the simulation. This 


practical exercise is what remains to be examined in this Appendix. 


4. COMPARISON OF JANUS SIMULATION RESULTS 
TO THE DETERMINISTIC LANCHESTER SQUARE 
LAW RESULTS 


Using the values for & and 6 derived from the Janus simulation data output, 


one can form the two Lanchester Square Law (F|F) Equations C.10 and C.11 below. 


d 

= = —(0.00217256)y(t) with zp = 113 (C.10) 
d 

rs = —(0.00072153)z(t) with yo =55 (C.11) 


Or, equivalently, using Equations C.5 and C.6 the force size equations can be written 


as 


where 


oS 


1 x a “a .* 
5 fn _ $55] exp v8) ai [u + (555 vp Vint (C.12) 


(55 = zu] exp(V)t 4 (ss = us) exp Vet (C.13) 


Qo] o 


a 0.00217256 


a 
| 


o~ 
I 


Q.00072153. 
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Calculating the force size values of x(t) and y(t) for values of t = (0,15, 30,--- , 210) 
minutes using Equations C.12 and C.13 above, the reader will note that the deter- 
ministic approximations depicted in Table XV on page 105 of this Appendix are quite 
representative of the Janus simulation results. (Note: The deterministic results for 
x(t) and y(t) are shown as integer results which were derived by rounding the contin- 
uous solutions to Equations C.12 and C.13). Employing the deterministic Lanchester 
equations to replicate the high resolution simulation results is called “replaying” the 
simulation. This is the role that the ATCAL Phase II algorithm plays as part of the 
aggregated models it supports. In this example, the deterministic Lanchester results 
estimate the casualties achieved during the Janus simulation with fairly high accuracy 
as depicted by Figures 11 and 12 on page 106. Note that the maximum absolute error 
of the approximations is no greater than two weapon systems for both forces. 

The small discrepancies can be attributed to several facts. First, the author ar- 
bitrarily decided to single out two opposing weapon systems from a highly integrated 
heterogeneous combat simulation to demonstrate the mathematics associated with the 
homogeneous form of Lanchester’s Square Law (F|F’) and the COMAN MLE attrition 
coefficient estimation process. The reader may have already gathered that the high 
resolution simulation presented in Chapter II infers the need for a heterogeneous-type 
attrition analysis construct; however, tractable formulations do not exist for het- 
erogeneous Lanchester-type equations. So, in presenting the simplistic closed form 
solutions associated with the homogeneous problem, the reader is able to gain an 
appreciation for the mathematics behind attrition computations. Second, examin- 
ing the casualty times listed in Table XIV, one can see that the exchanges of aimed 
fire between the REDFOR BMP-2 (with AT-5) and the BLUEFOR M2 (with TOW) 
were focused around several small battles. Each force shared losses during several 
distinct time periods, i.e., intervals where [t, — t,-)] are small and follow each other 
closely. The larger values of [t, — tz_] depict periods of inactivity during the battle. 
Perhaps distinct COMAN MLE attrition rate estimates @ and b could be derived for 
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several time sub-intervals representing their corresponding sub-battles. Lastly, the 
reader should recall that this research only utilizes one stochastic sample of the high 
resolution simulation. A more statistically robust analysis could have been made 
had the attrition rate parameters been estimated after considering the time series 
Markov (Poisson) attrition processes for several independent simulations. 

Even with the perceived shortcomings, the reader will agree that the sampled 
deterministic Lanchester replay of the high resolution simulation is quite good. Ad- 
ditionally, through the application and understanding of this simplified example, the 


reader should have a thorough appreciation for the expanse of this subject matter. 
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Casualty 


Time 


16.41 
16.56 
17.42 
18.27 
18.28 
37.26 
37.49 
38.10 
38.11 
46.32 
49.24 
99.29 
64.14 
83.12 
84.20 
84.21 
110.58 
ATE 3 
131.51 
142.28 
142.30 
143.03 
143.10 
149.26 
149.30 
149.48 
153.40 
173.36 
183.35 


[Ref. 4] 


Target 


No. 
Lost 


I 
] 
1 
I 
I 
| 
| 
1 
1 
] 
] 
I 
i 
] 
] 
] 
I 
] 
1 
1 
] 
1 
1 
1 
I 
1 
1 
1 
i! 


Firer 


k 


28 
2 


Time Interval 


Phe 


16:41 
0:15 
0:46 
0:45 
0:01 

18:58 
0:23 
0:21 
0:01 
Al 
Ry 2 

10:05 
4:45 

18:58 
1:08 
0:01 
Dok 
0:15 

20:38 

10:37 
0:02 
0:33 
0:07 
6:16 
0:04 
0:18 
S202 

19:56 
9:59 


(™mz-1) 


113 
EZ 
112 
111 
110 
110 
109 
108 
108 
107 
106 
106 
106 
106 
105 
105 
104 
103 
102 
101 
101 
100 
100 
100 

99 

ao 

98 

97 

96 


Table XIV. Janus Simulation Data: Coroner’s Report 
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No. of M2 No. of BMP-2 


(nx-1) 


55 
50 
54 
54 
54 
53 
53 
53 
52 
52 
52 
51 
50 
49 
49 
48 
48 
48 
48 
48 
47 
47 
46 
45 
45 
44 
44 
44 
44 


t 


Janus x(t) 


113 
I 
109 
108 
106 
105 
103 
101 
100 


Table XV. Comparison of Janus Attrition Data to Lanchester-Type Replay 


Lanchester x(t) 


113 
113 
110 
107 
106 
106 
104 
104 


Janus y(t) 


105 


Lanchester y(t) 
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